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We study the effect of critical fluctuations on the (B,T) phase diagram in extreme type-lf 
superconductors in zero and finite magnetic field. In zero magnetic field the critical fluctuations are 
transverse phase-fluctuations of the complex scalar Ginzurg-Landau order parameter, which when 
excited thermally will induce line-defects in the form of closed vortex loops into the system. The 
distribution function D{p) of vortex loops of perimeter p changes from an exponential function 
D{p) ~ p"" exp(— e(T)p/fcsT) to a power law distribution D{p) ~ p~" at the zero-field critical 
temperature T = Tc. We find that the long-wavelength vortex-line tension vanishes as e{T) ~ 
|r — TcP; 7 ~ 1.45, as r — > Tc- At T = Tc, an extreme type-II superconductor suffers an unbinding 
of large vortex loops of order the system size. When this happens, the connectivity of the thermally 
excited vortex-tangle of the system changes abruptly. When amplitude fluctuations are included, it 
is shown that they are far from being critical at the superconducting transition temperature Tc . The 
vortex-loop unbinding can therefore not be reparametrized in terms of critical amplitude fluctuations 
of the original local Ginzburg-Landau order parameter. The loss of phase-stiffness in the Ginzburg- 
Landau order parameter, the anomaly in specific heat, the loss of vortex- line tension, and the 
change in the connectivity of the vortex-tangle are all found at the same temperature, the critical 
temperature of the superconductor. At zero magnetic field, unbinding of vortex-loops of order the 
system size can be phrased in terms of a global [/(l)-symmetry breaking involving a local complex 
disorder field which is dual to the order parameter of the usual Ginzburg-Landau theory. There is 
one parameter in the theory that controls the width of the critical region, and for the parameters we 
have used, we show that a vortex-loop unbinding gives a correct picture of the zero-field transition 
even in the presence of amplitude fluctuations. A key result is the extraction of the anomalous 
dimension of the dual field directly from the statistics of the vortex-loop excitations of the Ginzburg- 
Landau theory in the phase-only approximation. In finite magnetic fields, the first order vortex-line 
lattice (VLL) melting transition is accompanied by a loss of longitudinal superfluid stiffness; this 
is true also for the isotropic case. A scaling analysis of the vortex lattice melting line is carried 
out, yielding two different scaling regimes for the vortex lattice melting line, namely a high-field 
scaling regime and a distinct low-field ?>DXY scaling regime. We also find indications of an abrupt 
change in the connectivity of the vortex-tangle in the vortex liquid along a line Tl{B), which at 
low enough fields appears to coincide with the VLL melting transition line within the resolution of 
our numerical calculations. We study the temperature at which this phenomenon takes place as a 
function of system size and shape. Our results show that this temperature decreases and appears to 
saturate with increasing system size, and is insensitive to aspect ratios of the systems on which the 
simulations are performed on, for large enough systems. A necessary, but not sufficient, condition 
for the vortex-line tension to vanish is a change in the connectivity of the vortex tangle in an extreme 
type-II superconductor. 
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I. INTRODUCTION 

Ten years after Abrikosov's classic prediction of a 
lattice of quaiitized vortices, the Abrikosov vortex-line 
lattice (VLL),EJ as the ground state of type-II super- 
conductors when the magnetic field is tuned beyond a 
lower critical valueH, Gerd Eilenberger suggested that the 
VLL caiild melt close to the critical temperature of the 
systemtl. The magnetic field versus temperature {B,T)- 
phase diagram of extreme type-II superconductors has 
for some time been under intense investigation both the- 
oretically and experimentally, following suggestions that 
the VLL could undergo a melting transition in regime 
of the {B, r)-phase diagram of the high-temperature aup 
perconductor that could be experimentally resolvedQ'B. 



This was soon confirmed by a more thorough theoretical 
analysia3 where it was shown that the VLL of the high- 
temperature superconductors was particularly suscepti- 
ble to thermal fluctuations due to the large anisotropy of 
these compounds. The anisotropy only affects the melt- 
ing line of the VLL when the pronounced nonlocal elas- 
tic properties of the VLL in strong type-II superconduc- 
tors, first discussed for the isotropic case in the pioneering 
works of BrandtLrB, are taken into accountQ. 

While it iJfOw appears well established both 
experimentallyg an d theoretically for three dimensional 
vortex systemst2ltll that the VLL in the clean limit of 
type-II superconductors melts in a first order phase- 
transition, much less consensus has been reached on how 
to describe the state which the VLL melts into, even 
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in the clean limit. Only very recently haajt been es- 
tablished, through numerical simulationst30 that the 
vortex-liquid is always incoherent, i.e phase-coherence is 
destroyed in all direction, including the direction of the 
induction, as soon as the VLL melts. Inside the vortex- 
liquid regime, there is therefore no transition from a 
disentangled to an entangled vortex-liqiud. For such a 
transition to occur inside the vortex-liquid, the longitu- 
dinal superfluid density would have to be non-zero above 
the melting temperatjuce- This however does not hap.- 
pen in the clean limitO'Ej, even in the isotropic caseliJ. 
Recently, questions have also been raised whether the 
vortex-line picture of the molten phase o£-tbe | A|br i kQ sov 
VLL is viable at all at low fields B < ITHlBiMiM. 

In terms of fundamental physics, extreme type-II su- 
perconductors are interesting due to their large fluctu- 
ation effects not commonly seen in condensed matter 
systems. This is ultimately due to the fact that they 
are strong coupling superconductors arising out of doped 
Mott-Hubbard insulators. The latter fact gives rise to 
the effect that the phase-stiffness of the superconducting 
order parameter is low, due to a low value of the super- 
fluid density ps 

where /S.9 is a twist in the superconducting order param- 
eter over the size of the system, and / is the free energy 
density. This particular and important aspect of doped 
Mott-Hubbard insulators has-beep quite strongly empha- 
sized already for some timeE2rEil, see also Ref. ^ The 
strong coupling effect gives rise to a large Tc, so that the 
Ginzburg-Landau parameter k ^ Td -Jpl is large. This 
also softens the vortex matter in these systems, garticu- 
larly when coupled with their strong layerednesaj. 

There is also a close connection between thermody- 
namic phase transitions |4ji these systems^ and phaser 
transitions in sunerfluidsE^, liquid crystalscj, crystalsEs, 
and cosmologycTEj. The close connection between these 
apparently different physical problems, is due to the sim- 
ilarity of the topological objects that appear in these 
problems. Particularly in the context of superfluid He^, 
the proposition that an unbinding of topological phase- 
fluctuations in the form of vortex- loops is the microscopic 
mechanism for the superfluid-normal ptpitei tpmsition, has 
been extensively studied in the pastEiTEJ and early 
attempts at formulating a field-theory of this in the con- 
text of charged superfluidsriji zero magnetic field has also 
appeared in the literatureEZl. Effective gauge-field theo- 
ries with an internal [/(l)-symmetry all have in common 
that they support line- deject in the form of vortex-loop 
excitations as stable topological objects. Understanding 
the role of such excitations on the (B, T) phase-diagram 
of type-II superconductors is an important problem in 
physics, and presumably will shed light on the related 
problems mentioned above as well. 

In conventional low-temperature superconductors, the 
temperature where Cooper pairs start to form, Tmf, is 



practically identical to the true superconducting transi- 
tion temperature Tg. The commonly applied Ginzburg- 
criterion provides a useful estimate for the width of the 
critical regions in systems with weak fluctuation effects, 
showing that the width of the critical region is of order 
\t\ - {Tc~T)/Tc - 10^6-10"''. A mean field description 
of the S-N phase transition is appropriate 'for all practi- 
cal purposes" . In high- Tc superconductors, this may no 
longer be the case. There appears to be mounting ex- 
perimental evidence that the width of-the, critical region 
is as large as a few Kelvin in YBC0lJl3, which would 
encompass the melting line of the flux-line lattice up to 
a field of order ITQ. 

In zero field the superconducting-normal phase 
transition.— js ,■ exclusively caused by a vortex loop 
unbindingl3'ollJa'C3. Below the critical temperature Tc 
vortex loops are confined to a typical perimeter Lqj ^md 
cause only local disturbances in the macroscopic super- 
conducting state. Recently, this has been demonstrated 
clearly, by correlating an abrupt change in vortex tangle 
connectivity, a loss of vortex-line tension, loss of super- 
fluid stiffness and specific heat anomaly precisely at the 
critical temperature of the superconductor, even for the 
isotropic caseCJO. At Tc , thermally induced vortex loops 
loose their effective line tension and therefore unbind. 

In this scenario, at low fields, thermally induced vor- 
tex loops could conceivably interact strongly with field 
induced flux lines. This interaction is ignored in mod- 
els using the line-only approximation, where the ther- 
mally induced vortex loops degrees of freedom are ne- 
glected, effectively being considered as irrelevant rela- 
tivistic correctioos in a corresponding 2D quantum bo- 
son system.B'Oc3. This model has met with consider- 
able success in describing parts of the flux-line lattice 
melting curve at intermediate to elevated fields, where 
its position the phase-diagram as well as its depen- 
dence oruanisotropy was explained using the Lindemann- 
criterioiiaE3c3. The melting line for fields of more than a 
few Tesla is little affectedly the vortex-loop unbinding, 
as pointed out recentlyCj'CJ. 

However, the question arises whether this is a tenable 
conclusion for low fields as well. The fact that the zero- 
field transition can be characterized precisely by a loss of 
line-tension of thermally induced flux lines, implies that 
there is a sharp change in the distribution function D{p) 
for vortex loops of a given perimeter p. It changes from 
an exponential form D{p) ~ P~r^ G^p{~~£{T)p/kBT) to 
a power law D{p) ~ at tJI^. This implies thCje^is.- 

tence of a diverging length scale LoiT) = kBT/e{T)MM. 
Given this fact, it raises the question of whether the 
critical fluctuations can affect the melting line in a size- 
able field-temperature regime, rendering the vortex lines 
tension-less. The vortex line tension is analagous to the 
mass of the bosons in a 2D non-relativistic boson-analogy 
of the vortex system. If the vortex-line tension were to 
vanish, it would mean that the boson-mass wouild vanish 
in the corresponding analogy. There is no non-relativistic 
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limit of any mass-less theory. The conclusion would be 
that any 2D boson- model which is non-relativistic, is in- 
applicable in the part of the phase-diagram where the 
vortex-line tension vanishes. We reemphasize that at ele- 
vated fields, where the first order flux- line lattice melting 
line splits off from the new transition line proposed here, 
the Lindcmann-criterion of flux-line lattice meltings is 
ex;^^ted to correctly locate the position of the melting 

The outline of this paper is as follows. In Section II, we 
introduce the Ginzburg-Landau model studied in this pa- 
per, and various approximations and reformulations of it, 
as well as their inter-relations. In Section III we present 
the ideas underlying the simulations that are presented 
in this paper, and introduce and discuss the quantities 
we study. In Section IV, we present results of the simu- 
lations in zero magnetic field. In particular, we present 
results which demonstrate that the zero-field transition in 
an extreme type-II superconductor is driven by a prolif- 
eration of unbound vortex loops, which therefore consti- 
tute the critical fluctuations of this system. In Section V, 
finite-field results are given. Summary and conclusions 
are presented in Section VI, and in this section we also 
list point by point the new results obtained in this paper. 



II. MODELS 

In this section we define the models considered in this 
paper: 1) the continuum Ginzburg-Landau model, 2) the 
lattice Ginzburg-Landau model in a frozen gauge approx- 
imation, and 3) the uniformly frustrated 3D XY model. 
We also discuss the approximations involved and the va- 
lidity of the models. 



/i-direction, $o = h/2e is the flux quantum, and /iq is 
the vacuum permeability. In Eq. 0, the gauge field A is 
related to the local magnetic induction, 6(r) = V x A(r). 
Finally, the GL parameter g is assumed to be temper- 
ature independent, while a = o(T) changes sign at a 
mean field critical temperature Tmf{B), where Cooper 
pairs start to form. B is the spatial average of the 
magnetic induction. The critical temperature Tc where 
phase-coherence develops, is always smaller than Tmf', 
the existence a finite Copper-pair density does not imply 
that the system is in a superconducting state. 

Later on, we shall recast the Ginzburg-Landau the- 
ory in a quite different form that also exhibts a U{1)- 
symmetry, but where the field conjugate to the relevant 
phase is the number operator for the topological excita- 
tions destroying the order of the Ginzburg-Landau theory 
itself. Although this may seem like an unnecessary com- 
plication, it has the advantage of facilitating a detailed 
discussion of the vortex-liquid phase of the GL-theory 
in terms of the ordering of some local field, namely the 
complex scalar field (f){r) to be introduced and discussed 
in Section IIF. This is not possible using the Ginzburg- 
Landay function, i/jir), sin ce . < ipi'"^) > is always zero 
in the vortex liquid phasetBEj. In the zero-field low- 
temperature ordered phase, the system spontaneously 
chooses a preferred phase angle Q, and explicitly breaks 
the U{1) symmetry. The vortex-sector of the GL-theory 
also exhibits a ?7(l)-symmetry breaking, but where U{1)- 
symmetry is broken in the high-temperature phase, and 
restored in the low-temperature phase. 

Eq. |l| has two intrinsic length scales, the mean-field 
coherence length 



2m^\a{T)y 
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A. Ginzburg-Landau model 

Our startitig point is the continuum Ginzburg-Landau 
(GL) modeH. In quantum field theory, the GL model is 
also referred as the scalar QED or the J7(l)-f Higgs model 
or the Abelian Higgs model. The effective Hamiltonipj 
for the GL model in an anisotropic system is given bjO 



Hgl = I d?r 



a{T) + 



27r 



2m„ 



^ $0 



(1) 



Here, = |'0(r)|e*^^''' is a complex order field rep- 

resenting the superconducting condensate. In supercon- 
ductors, the amplitude |V'(r)P should, be interpreted as 
the local Cooper-pair density. Furthermore, m^^ is the ef- 
fective mass for one Cooper pair when moving along the 



and the magnetic penetration depth 
^ 4/ioe2|a(r)|- 



(3) 



is the characteristic length of the variation of |^/'(r)| 
along the /^-direction, and is the characteristic length 
of the variation of the current flowing along the /x- 
direction. 

In order to carry out Monte Carlo simulations of the 
GL model, the model is discretized by replacing the co- 
variant derivative in the continuum GL Hamiltonian, Eq. 
|l|, with a covariant lattice derivative, 

27r 

a^ \ / 

The resulting model is a version of the Lawrence-Doniach 
modeH with all three directions discretized. The effec- 
tive Hamiltonian for the lattice GL model is given by. 
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E 



«lV'P + f|^|' 



fi—x,y,z 



V'(r + /i)e-**^'^''^^W -^(r) 



fi—x,y,z 



(5) 



Here, and /i is a lattice constant and a unit vector 
along the //-axis, respectively. Furthermore, the lattice 
derivative is defined as 

A^V(i") = + A) - V'(r)- 

Taking the continuum limit — > 0, the effective Hamil- 
tonian for the lattice GL model (Eq. ||) reduces cor- 
rectly to the GL effective Hamiltonian in the continuum 
(Eq. H). As defined in Eq. ||, the lattice GL model does 
not contain vortices. To reintroduce the vortices in the 
model, we must compactify the gauge-theory by requir- 
ing that the gauge invariant phase differences satisfjCj, 

[0(x + A) - e(x) - ^aM^)] e [-^, ^ > ■ (6) 

Whenever this constraint is used to bring the gauge in- 
variant phase differences back to its primary interval, we 
automatical! introduce a unit closed vortex loop, and the 
net vorticity of the system is guaranteed to be conserved 
at every stage of the Monte-Carlo simulation. From the 
renormalization group point of view the continuum GL 
model and theJattice GL model belong to the same uni- 
versality clas£3. We therefore expect the lattice GL 
model and the continuum GL model to give, qualita- 
tively, the same results. 



B. Lattice Ginzburg-Landau model in a frozen gauge 
approximation 

In extreme type-II superconductors, the zero temper- 
ature mean-field penetration depth is much greater than 
the zero temperature coherence length, A^(r = 0) ^ 
£,fi(T = 0). Thus, fluctuations of the gauge field repre- 
sented by the last term in Eq. ^ around the extremal 
field configuration are strongly suppressed and can there- 
fore be neglected. The effective Hamiltonian for the 
frozen gauge (FG) model can be written as 



rr l"(0)P 
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^—x,y,z 



^ |-|V'(r + A)IIV''(r)| [2-2cos(A^0 -A,)] 



(7) 



Here, we have defined a dimensionless order field and 
vector potential 



j(0)| 



27r 







The natural energy scale along the /^-direction is, 

t2 
-,2 



9 al 



Assuming a uniaxial anisotropy along the z-axis, the nat- 
ural energy scale for the FG model is 



— Jx 



2|a(0)P 



(8) 



Here, we have put our coordinates {x,y, z)-axis parallel 
to the crystals (a, b, c)-axis. Furthermore, = = Cafe 
and = is the coherence length in the CuO-planes 
and along the crystal's c-axis, respectively. Furthermore, 
Aa; = Aj, = Xab and A^ = Ac is the penetration depth in the 
CuO-planes and along the crystals c-axis, respectively. 
In Eq.^, d is the distance between two CuO supercon- 
ducting planes in adjacent unit cells. The energy scale 
Jo is rpugbly the energy scale of exciting a unit vortex 
IoopBHH. 

The ratio between the energy scales and serves 
as an anisotropy parameter. 



r = 



^c.O'x ^ab^x 



(9) 



In this model, the lattice constant should be defined 



a,, = max(d^, Co^^). 

Here, d^j is an intrinsic length along the /^-direction in 
the underlying superconductor to be modeled. Examples 
of such intrinsic length are the distance between CuO- 
planes in adjacent unit cells, the (a,b)-dimension of the 
unit cell. To be consistent, the constant Co should be 
larger than ^ 4. This requirement a^/^fj, > 4 ensups 
that the amplitude of the order field does not overlapEj. 
Such overlap will give rise to a domain wall term (V|';/'|), 
which is absent in the lattice GL model. 

Within the frozen gauge approximation, the gauge field 
serves only as a constraint in fixing the value of the uni- 
form induction. In terms of magnetic induction this ap- 
proximation is valid when B 3> Bci{T), where the field 
distribution from individual flux lines overlap strongly, 
giving uniform induction. Note that Bci(T) also van- 
ishes when the temperature approaches Tc- In zero field, 
this approximation is valid for all temperatures except 
an inaccessible temperature region around rj£3. 

In our simulations on the FG model, we allow for both 
phase- and amplitude fluctuations of the superconducting 
order parameter tp{r) = |V'(r)| exp[j6'(r)]. Details of the 
Monte-Carlo procedure for this case will be given below. 
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C. Uniformly frustrated 3D XY model 

The uniformly frustrated 3D XY model was first used 
as a phenomenological [iHodel for extreme tepe-II super- 
conductors by Li et alO and Hetzel et al.til. To obtain 
the uniformly frustrated 3D XY (3DXY) model from the 
FG model, we freeze the amplitude of the complex order 
field in Eq. \tp'\ = 1. This is the London approxima- 
tion. The resulting effective Hamiltonian for the 3DXY 
model is given by 



H 



XY 



2|a(0)|2 , , 
axttya^ cos (A^fc* - A^j 



(10) 



The lattice constants in the 3DXY model should be de- 
fined as 

= max(d^,^^). 

Assuming uniaxial anisotropy, the energy scales and the 
anisotropy parameter of the 3DXY model are the same 
as for the FG model, Eqs. | §. Note that both the FG 
model and the 3DXY model contain precisely the same 
topological objects, i.e. vortex loops and vortex lines, as 
for the GL model. The local gauge symmetry in the GL 
model is however reduced to a global f7(l)-symmetry in 
the FG model and the 3DXY model. 



D. Villain-approximation and vortex representation 

To further corroborate interpretations of the results 
from our Monte-Carlo simulations using the uniformly 
frustrated 3DXY-mode\, to be detailed in the next sec- 
tion, it is useful to provide an alternative, but entirely 
equivalent formulation of the GL-theory. This formula- 
tion replaces a description in terms of the GL-function 
-0 by vortex-degrees of freedom, where the interaction 
between vortex-segments is mediated by a gauge-field, 
which we denote by h. This gauge-field is not the elec- 
tromagnetic vector potential A, but will couple to it. 
The resulting structure of the theory makes it possible, 
in three dimensions, and three dimensions only, to refor- 
mulate the vortex-content of the GL-theory as a theory 
of a complex matter field </) coupled to the gauge-fields h 
and A. Although this may seem as an unnecessary de- 
tour, the great advantage of this approach, is that certain 
vortex- correlators, notably our quantity Ol to be defined 
below, can be directly related to a C/(l)-symmetry of the 
(/)-theory. 

To proceed with this, we introduce the well-known 
Villain-approximation to the 3DXY model. The Villain 
approximation consists of replacing the cosine potential 
in the uniformly frustrated 3-DA"Y-model by a Gaussian 
27r-periodic potential. In this way the longitudinal spin- 
wave excitations of the 0-field decouple from transverse 
vortex-excitations of the theory. This decoupling does 



not alter the critical behavior of the system. The parti- 
tion function for this theory reads, after a rescaling of 
the vector potential and charge 



Zv^li J De{r) J DA{r) ^1 

r n(r) — — oo 

S = - Etf - 2eA - 2nnr + i(V x A)^], 



where n(r) is an integer-valued field. The kinetic term 
is linearized by a Hubbard-Stratonovich decoupling, in- 
troducing an auxiliary velocity-field v(r) and using the 
identity 



2P 



.(r).u(r)]). 



This is now inserted back into the original partition func- 
tion, using u = V6' — 2eA — 27rn(r). The sum over the 
integers n(r) may be carried out, yielding the constraint 
that v(r) is integer valued, say v(r) = l(r). The next 
step is to integrate out the 0(r)-variable, which yields 
the constraint V • v(r) = 0, which is solved by introduc- 
ing an integer valued field h(r) such that l(r) = V x h(r). 
In order to be able to treat h as a continuous variable, 
we introduce a new integer-valued field m and apply the 
Poisson-summation formula 



E 



m— — oo 



E '^(i-h)' 



l=-c 



Note that this procedure does not involve any approxi- 
mations. Finally, we write the partition function for the 
GL-theory in phase-only and Villain-approximations as 



/oo 
DhDA E e^-I™''^^^!, 



(11) 



5cff - - E 

r 

1 



h + — (V X h)2 + 2ie(V X h) • A 
Zp 



-(VxA) 



where the following constraints apply in the functional 
integral: V-A = V- h = V- m = 0. The effective action, 
Eq. |ll|, is invariant under 



h 
A 



h + Vcjh 
A + VwA. 



The field h is readily interpreted as a fictitious gauge-field 
that mediates an interaction between vortex-segments m. 
This is easily seen by integrating out the h-ficld. 



5 



E. Dual (disorder field) representation 

Whenever a field theory sustains topological defects, 
it is often useful to formulate a new field-theory of the 
topologial excitations of the original theory per se, and 
this forms a dual description of the original theory. We 
will do this for the present problem also, following Ref. 

This means that the vortex-content Eq. |l^ of 
the Ginzburg-Landau theory Eq. |^ in the phase-only 
Villain-approximation is cast into the form of a local 
field theory involving a complex scalar mass-field describ- 
ing local vortex-fluctuations, coupled to a dual gauge- 
field that mediates an interaction between the vortex- 
segments. The resulting theory will exhibit explicitly a 
C/(l)-symmetry, and as always in such cases, the ques- 
tion to be asked is under what circumstances, if any, the 
symmetry will be spontaneously brokencZI. 

The purpose of this reformulation is to provide a point 
of contact between on the one hand a quantity to be 
introduced in Section IIIA and studied in Sections IVB 
and VB, probing vortex-tangle connectivity and denoted 
Ol, and on the other hand thermodynamics. The key 
point is that in zero magnetic field, the two-point cor- 
relation function of the complex scalar mass-field 0(r) 
of the dual theory, is precisely the probability of finding 
a connected vortex path between the two points of the 
correlation function, |*egardless of by which path the two 
points are connectedc3. Long-range order in G(x, y) im- 
plies a broken J7(l)-symmetry. Equivalently, long-range 
vortex connectivity in zero magnetic field implies a bro- 
ken U {l)-symmetry, which is "hidden" at the level of Eq. 
pr] , but is brought out when Eq. |ri| is reformulated to 
the dual form, to be described below. 

In three dimensions, and three dimensions only, a 
vortex-loop system interacting with a long-range Biot- 
Savart interaction and steric repulsion, may in the con- 
tinuum limit be written as a gauge .theory of a local com- 
plex matter field (f> coupled to hcjOEa. We may extend 
the results of this work including fluctuations in A in a 
finite magnetic fleld, in which case the vortex-content of 
the Villain-approximation to the GL-theory corresponds 
precisely to an action of the following form 

Z = YlJ D(t>{Y) D(t>*{T) Dh{v)DA{Y)e^^''^'>'''>''^^^^^ 



5, 



off 



E 



"Vl' + ^l0l' + Jl(--e'h)0|^ 
2 2, I 



the present (7(l)-symmetric theory is the dual description 
of the Ising- model, where a local fleld is introduced to de- 
scribe the local fluctuations in the topological defects of 
that model, which are domain walls separating different 
spin-ordered regions of the magnet. The effective action, 
Eq. |l^, is invariant under the set of transformations 

(j) 0exp(ia;h) 



^(V X h)2 + 2i e (V X h) • A + i(V x A)^ 



, (12) 



where the coefficients {a',e',g') appearing in the theq 
are given in terms of the parameters entering Eq. pE 
For our discussion, their precise values are of no impor- 
tance. The interpretation of the 0-field is that it is a local 
field describing local fluctuations in the topological exci- 
tations of the GL-theory, namely line-defects in the form 
of vortex-lines. An analogue of this dual description of 



h 

A 



h + 4 Vwh 

e' 

A + VcjA- 



(13) 



By rewriting the theory in Eq. in this way one ob- 
serves that it explicitly exhibits a f7(l)-symmetry. Note 
that this relies entirely on the possibility of reformulat- 
ing the interacting loop-gas, including Biot-Savart inter- 
actions, in terms of a complex matter field (j) coupled to 
a gauge-field h. 

Consider now for the moment the case of zero magnetic 
field. The main advantage of the above formulation is 
that the probability of finding a connected path of vortex 
segments, starting at x and ending at y, G'(x, y), is-given 
by the two-point correlation function of the ^-fielda 



G(x,y) r(x) 0(y)> 



(14) 



A vortex-loop unbinding will lead to a finite G(x, y) when 
|x — y| — > oo, because infinitely large loops will connect 
opposite sides of the vortex system. On the other hand, if 
lim|x_y|^oo G'(x, y) 7^ 0, this implies that < (/)*(x) >^ 0, 
corresponding to a broken J7(l)-symmetry. Therefore, 
the dual field 0(r) is an order parameter of a vortex-loop 
unbinding transition. The broken t/(l)-symmetry is asso- 
ciated with the loss of number conservation of connected 
vortex-paths threading the system in any direction (in- 
cluding direction perpendicular to an applied magnetic 
field, if that is present). This limit of the two-point cor- 
relation function is closely related to the quantity Ol we 
introduce in Section IIIA, which probes the connectivity 
of the vortex tangle in an extreme type-II superconduc- 
tor. The above connection makes it at the very least 
plausible that an abrupt change in this connectivity, as 
probed by the change in Ol, is associated with breaking 
a f7(l)-symmetry of the ?;orte2:-sector of the GL-theory, 
equivalently an onset < 4>* > 01 < 4> >. Since this only 
happens above a critical temperature, we may view as 
a disorder-iield, in contrast to the order-parameter field 
i/j of the original GL-theory. We will make explicit use of 
this connection in Section IV. 

In zero magnetic field, we will show that the loss of su- 
perfluid density, specific heat anomaly, change in vortex- 
loop distribution, loss of long-wavelength vortex-line ten- 
sion, and abrupt change in vortex tangle connectivity 
abruptly occurs precisely at the same temperature both 
for the 3DA"F-model, also when amplitude fluctuations 
are included. 

At flnite magnetic fields, the situation is complicated 
by the fact that the vortex system is always connected 
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across the system in at least one direction, namely the 
field direction, at all temperatures. One may how- 
ever still extract information of the type encoded in 
< (j>* {r)(j){r) > at zero field by performing a singu- 
lar gauge-transformation of the type used in Ref. ^ 18, 
which roughly speaking amounts to subtracting out the 
field-induced vortices and studying the remaining loop- 
gas, which has a field-theory description very much like 
the zero-field version of Eq. |l^. The obvious advantage 
of this is that one removes the asymmetry of the system 
imposed by the magnetic field. A twopoint correlator 
of this theory then probes the connectivity of non-field 
induced vortex paths across the system, which in turns 
probes the possibility of having a broken t/(l)-symmetry 
and hence an onset of the order parameter < (f){r) >^ 0. 

We will perform a numerical analogue of this in 
our simulations, namely we will probe the connectiv- 
ity of the vortex tangle of the superconductor in direc- 
tions perpendicular to the magnetic field. Ideally, what 
one should do is to generate phase-configurations (and 
vortex-configurations) of the extreme type-II supercon- 
ductor, subtract out from each configuration a number 
of vortex paths that connects the system along the field 
direction precisely corresponding to the number of field 
induced vortices in the system, which is a fixed number 
in a canonical ensemble usually studied for this problem. 
Out of the remaining vortex tangle one may then try to 
find vortex paths connecting opposite sides of the system. 
Numerically this procedure is entirely prohibitive and we 
therefore opt for the algorithm of calculating Ol, to be 
described in detail in Section IIIA. 

We stress that the procedure of computing Ol de- 
scribed in Section IIIA unquestionably probes the con- 
nectivity of a vortex-tangle across the system, not asso- 
ciated with magnetic field, precisely as in the zero-field 
case. The objective is to probe the breaking of a U{1)- 
symmetry associated with the proliferation of unbound 
vortex-loops in the system, as pointed out in Ref. 
This will be shown to be precisely borne out in zero 
magnetic field. In finite magnetic field we also obtain 
a weak specific heat anomaly at the temperature where 
Ol changes abruptly, as the system size is increased. 



III. DEFINITIONS, SIMULATION PROCEDURE, 
AND MODEL PARAMETERS 



In this section, we define the physical quantities con- 
sidered, describe our Monte Carlo procedure, and present 
the values of the model parameters used in the simula- 
tions. 



1. Specific heat C 

To calculate the specific heat per site C, we use the 
fiuctuation formula. 



C 1< > - < >2 



kB V 



(15) 



and is 



Here, the dimensionless volume V — L^LyL^ 
the dimensionless linear dimension of the system along 
the /i-direction. is measured in units of the lattice 
constant a^. As a check of consistency, we also calculate 
the specific heat per sit' 
of the internal energy 



the numerical derivation 



l_dU_ 
Vdf' 



(16) 



Note that for the FG model, where the effective Hamilto- 
nian depends explicitly on the temperature, there strictly 
speaking is an |additional term in the expression for the 
internal energyt3. 



U 



d\nZ 

d(3 



--<H> 



For the 3DXY model, if the simulations are properly 
done, C = Cu- For the FG model, where the effec- 
tive Hamiltonian explicitly depends on the temperature, 
C ^ Cu- Note however that < dH/dT > arises out of 
introducing a terappjfttfii'e dependence of the coefficients 
of the GL-theoryEllE3il3. The temperature dependence of 
these coefficients always has a temperature dependence 
set by the mean-field critical temperature. Thus, close to 
the true Tc these corrections, arising from integrating out 
the fermions of the iuiderlying microscopic description of 
the superconductorEil, are always negligible compared to 
the contribution coming from the true critical fluctua- 
tions of the order parameter, i.e. the transverse phase- 
fluctuations. In terms of Eq. ||the term — T < dH/dT > 
orginates from the temperature dependence of the ampli- 
tude of the order-parameter. Were this to actually van- 
ish dXT — Tc, substantial corrections to the specific heat 
and entropy would result. As we will see later, at T = Tc, 
the ensemble average of the amplitude of the order pa- 
rameter is far from renormalized to zero by vortex-loop 
fluctuations. Hence, at the critical point the correction 
term — T < dH/dT > in the .iiiternal energjfr-iwith its 
resulting corrections to entropyE^, is negligibleuJ. There 
is now ample experimental evidence that critical fluctu- 
ations are indeed important over a sizeable temperature 
window in || tho . hiffhr Tc cuprates of order several Kelvin 
below TcHEjC3Se3 a window which is consistent with a 
coherence length of order 10 A, about two orders of mag- 
nitude shorter than in conventional superconductors. 



A. Definitions 



2. Local Cooper-pair density < > 

As a probe for the local Cooper-pair density, we calcu- 
late 
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<|V''P>^-^<|V/(r)|^> 



(17) 



We see in Eq. that < > involves both thermal 

and space average. Recall that V' = V/\/|a(0)|/5- At 
the mean field level, we expect < > to develop an 
expectation value below the mean field critical tempera- 
ture TMi.(B). 



Jn 



ksTV 



For the FG case, 



3. Superfluid condensate density \ <tp' > |^ 

As a probe for the local condensate density (density 
of Cooper pairs participating in the superconducting con- 
densate), we calculate 



<^'>P^l5]|<V/(r)>|^ 



(18) 



Note the diff'erence between < > and | < tp' > p. 

The former describes local Cooper-pair density, while the 
latter describes what is commonly known as the conden- 
sate density in ^i7e-physics. In zero field, we expect 
I < t/y > p to develop an expectation value below the 
critical temperature Tc. 



^ = y ( E + A)l cos[A^0 - A,] 



Ju 



ksTV 

^ ^|^'(r)||V''(r + A)|sin[A^0-^, 

Note the difference between | < -!/;'> p and T^, they are 
not identical. The former quantity probes the superfluid 
condensate density, which is a locally defined quantity, 
while the latter quantity probes a global phase coher- 
ence along a given direction fj.. Since < V'' > is the order 
parameter of the Ginzburg-Landau theory, close to the 
critical point we have 



^. Distribution of the order field phase angle 

To probe the distribution of the phase angle in -0' (r) 
{r)\e^^'^^\ we define the distribution function 



(19) 



I < V'' > P ~ |t| 



2/3 



(20) 



where t — {T ~ Tc)/Tc. On the other hand, cx ps^, 
where is the superfluid density in the /i-direction. 
Using the Josephson scaling relation ps/j ~ 
|^|i/(ci-2)|^ along with the scaling laws 7 = i'{2 — rjjiJ 
and 2/3 = 2 - a - iB, we find 



Here, Sij is the Kronecker delta function. In the simula- 
tions, we have chosen to work with a discrete set of phase 
angles, 9', 9{r) = 2TTn/Ng. Here, n e [0, Ng] is an integer, 
and Ng is the number of allowed discrete phase angles. 
In our experience, the simulation results do not depend 
on Ng, provided Ng ~ 16. In zero field, when the phase is 
disordered, we expect Dg {9) to be uniformly distributed, 
D{9) = 1/Ng. In the ordered phase, we expect Dg{9) to 
show a peak around a preferred phase angle. 



5. Helicity modulus 



To probe the global superconducting phase coherence 
across the system, we consider the helicity modulus 
T^, defined as the second derivate of the free energy 
with respeGt|rtc) an infinitesimal phase twist in the p- 
directiorCJ-E^llZl. Finite means that the system can 
carry a supercurrent along the /i-direction. Within the 
3DXY-model, the helicity modulus along the /i-direction 
becomes. 



12/3-7, 1. 



(21) 



Here, d is the dimensionality of the system, v is the cor- 
relation length exponent of the system, (3 is the order pa- 
rameter exponent, 7 is the order parameter susceptibility 
exponent, and rj is the anomalous dimension of the order 
parameter two-point correlation function at the critical 
point. Therefore, although | < t/"' > p and are in 
principle different, they may appear to be very close if 
the anomalous dimension 77 of the ■(/'-field is small, as in- 
deed is ijbe case for the Ginzburg-Landau theory, where 
Tj ~ 0.04E3. Note that for ?7 > 0, the curve for should 
bend more sharply towards zero at the critical point than 
I < -0' > We will explicitly show by direct calcula- 
tions within the Ginzburg-Landau theory that | <ip' > p 
is very close to both in zero field and finite magnetic 
field. In zero magnetic field this is precisely what one 
would expect based on the above, when 77 << llI3. For 
the special case of d = 3, we have 2(3 ~ rjv = v < 2/3. 
To high precision, we have for the iDXY -model, that 
u = 0.673 and 77 = 0.038Ea 
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6. Vortex loop distribution D{jp) 

To probe the typical perimeter Lo{T) and the effec- 
tive long-wavelength vortex-line tension s{T) (not to be 
confused with the flux-line tension, which is always zero 
when gauge-fluctuations are completely suppressed due 
to the absence of tubes of confined magnetic flux), we 
define a vortex-loop distribution function D(p), which 
measures the ensemble-averaged numhct,|Of vortex loop 
in the system having a perimeter pOOalSa. In order to 
compare results from different system sizes, we normalize 
D(p) with respect to the system size. 

We search for a vortex loop using the following proce- 
dure. Given a vortex configuration, we start with a ran- 
domly chosen unit cell with vortex segments penetrat- 
ing its plaquettes. We follow the directed vortex path 
and record the trace. When the directed vortex path 
encounters a unit cell containing more than one outgo- 
ing direction, we choose the outgoing direction randomly. 
When the vortex path encounters a previously visited 
unit cell, i.e. when it crosses its own trace, we have a 
closed vortex loop, its perimeter being p. We now delete 
the vortex loop from the vortex configuration, to prevent 
double counting, and continue the search. The search is 
continued until all vortex segments are deleted from the 
system. 

Using a 3D non-interaeting boson analogy to the vortex 
system, it can be shownO that-the distribution-function 
D{p) can be fitted to the formO 



D{p) = A exp[- 



eiT)p 



(22) 



Here, A is a temperature independent coristant, and the 
exponent a w 5/2 to a first approximationZ^. When e(T) 
is finite, there exists a typical length scale Lq — ksT/e 
for the thermally excited vortex loops. The probability of 
finding vortex loops with much larger perimeter than 
is exponentially suppressed, according to Eq. When 
e = 0, D{p) decays algebraically, and the length scale of 
the problem Lq = kBT/e{T), has diverged. As a conse- 
quence, configurational entropy associated with topolog- 
ical phase-fiuctuations is gained without penalty in free 
energy. In zero field, there is only one critical point, and 
in this case Lq must be some power of the superconduct- 
ing coherence length £,{T). 



7. Probe of vortex- connectivity, Ol 

For probing the connectivity of a vortex tangle in a 
type-II superconductor, in zero as well as finite magnetic 
field, we introduce a quantity Ol, defined in zero mag- 
netic field as the probability of finding a vortex configu- 
ration that can have at least one connected vortex path 
threading the entire system in any direction. In the pres- 
ence of a finite magnetic field, Ol is defined as the prob- 
ability of finding a similarly connected vortex path in a 



direction transverse to the field direction, without using 
the PBC along the field direction. In zero field, we use 
the same procedure as in finite-field, namely searching for 
connected vortex paths perpendicular to the z-direction, 
although in this case we could just as well have used any 
direction. Note that Ol is very dijfereni^^Mm the wind- 
ing number W in the 2D boson analogy lljo. There, W 
is proportional to the number of vortex paths percolating 
the system transverse to the field direction. However, in 
the calculation of W, the PBC along the field direction 
is used many times. 

In an attempt to probe "vortex-percolation" , a slightly 
different quantity than Ol has been considered in 
the caatext of high-temperature superconductors by 
othersL3LZl. A crucial difference between our work and 
that of Ref. |7^, is that Ref. ^ allows periodic bound- 
ary conditions along the field direction to be used several 
times before the vortex path winds once around the x- 
or y-axis, as is easily seen from Fig. 2b of Ref. |7^. This 
ultimately is the same as computing the winding number 
of the 2D non-relativistic boson-analogy of the vortex- 
systemQ, as recently done in careful Monte-Carlo simu- 
lations in Ref. ^ It also explains why the authors of 
Ref. 1?^ get longitudinal dissipation at the onset of what 
they denote "vortex-percolation", which is nothing but 
the temperature at which the winding number becomes 
finite. 

This is entirely consistent with a number of 
other -MoiitepCaiJo simulation results on the 3DXY- 
modeO'OlijCJ'Eil which all show the loss of longitudi- 
nal phase-coherence and onset of longitudinal dissipation 
precisely at the vortex lattice melting transition. This is 
measured simply by the helicity modulus T^, which is 
quite different from Ol- To the contrary, in our calcula- 
tion of Ol, we do not allow for the use of periodic bound- 
ary conditions in the z-direction to measure vortex-tangle 
connectivity in the x- or y-directions, in other words the 
"percolating" configurations of Fig. 2b of Ref. ^ are not 
counted when computing Ol- 

We have 



Or = 



total 



(23) 



Ntotai is the total number of independent vortex config- 
urations provided by the Monte Carlo simulation. Fur- 
thermore, Nc is the number of vortex configurations con- 
taining at least one directed vortex path that traverses 
the entire system perpendicular to the direction, without 
using the PBC along the field direction. For convenience, 
we treat the zero field case as the limit VmiB~,o keeping 
the "field direction" intact. 

We search for the possibility of finding a vortex path 
such as described above by using the following procedure. 
Assume that the magnetic induction points along the z- 
axis. We follow all paths of directed vortex segments 
starting from all four boundary surfaces with surface nor- 
mal perpendicular to z, and check whether at least one 
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of these vortex paths percolates the system and reaches 
the opposite surface, without applying the PBC in the z- 
direction. Note that when crossing vortex segments are 
encountered, the procedure is to attempt to continue in 
a direction that will bring the path closer to the opposite 
side of the system, rather than randomly resolving the in- 
tersection. Ol is therefore a necessary, but not sufficient 
condition for finding an actual vortex-path crossing the 
system. However, in zero field this procedure does not 
make a difference to that of resolving the intersections 
randomly. This is demonstrated by the correlations of 
the change in and D{p), to be detailed in the next 
section. 

If a vortex path is actually found crossing the system 
in any direction in zero field, or without using PBC in 
the field-direction when a field is present, one may safely 
conclude that the vortex-line tension has vanished. If it 
were finite, it would not be possible to find such a path 
at all, either because all vortex-lines form closed confined 
loops in zero field, or because the vortex-line fluctuations 
along the field direction would be diffusive in finite field. 
In zero-field, this is clear by the above mentioned cor- 
relation between the change in Ol and D{p), cf. the 
results of the next section. In this paper, we also investi- 
gate this in detail for the finite-field case, by considering 
the position of the lowest temperature Tl where we have 
Ol = 1 both as a function of system size and aspect 
ratio Lx/Lz = Ly/L^. If vortex-line physics remains in- 
tact, Tl should move monotonically up with system size, 
and should scale with L^/Lz- Instead, we will find that 
Tl moves down slightly, and saturates with increasing 
system size at fixed aspect ratio. In addition, we find 
that Tl is virtually independent of aspect ratio for large 
enough systems. 

This contradicts expectations based on a lines-only ap- 
proximation to the vortex-liquid. It demonstrates that 
the connectivity of the vortex-tangle undergoes a funda- 
mental change inside the vortex-liquid. The above men- 
tioned finite size scaling analysis, suggests to us that this 
geometric transition is a property that survives in the 
thermodynamic limit. The issue is whether the change in 
connectivity has anything to do with a thermodynamic 
phase-transition. This will be investigated in detail for 
zero magnetic field in Section IIIA, and for finite mag- 
netic field in Section IIIB. In particular, we look for a 
specific heat anomaly scaling up with system size, at 
the putative transition point Tl- This will reveal if the 
change in the geometric properties of the vortex liquid is 
indeed associated with singular thermodynamics. In any 
case, once the geometric transition has taken place, it is 
no longer possible to model the vortex-liquid regime in 
terms of field-induced fiux lines only, with merely renor- 
malized interactions between them. 

In the VLL phase Ol = 0, since the field induced flux 
lines are well defined and do not "touch" each others, 
and the thermally excited vortex loons, are confined to 
sizes smaller than the magnetic lengthlill. Ol = 1 in the 
normal phase above the crossover region where the rem- 



nant of the zero field vortex loop blowout takes place. 
Needless to say, it is a matter of interest to investigate 
precisely where Ol changes value from zero to one. 

Note that Ol itself is not a genuine thermodynamic 
order parameter, althougk it may be said to probe an 
order-disorder transitionCJ. However, by the transcrip- 
tion of the vortex-content of the Ginzburg-Landau the- 
ory to the form Eq. ^ in Section IIF, it is brought out 
that probing the vortex-tangle connectivity by consider- 
ing Ol is closely connected to probing the two-point cor- 
relator of a local complex field <j){r), the dual field of the 
local vorticity-field m^(r) of the Villain- approximation 
and London-approximation to the Ginzburg-Landau the- 
ory, Eqs. The two-point correlator < 0*(r)(/)(r') > 
is ultimately the probe of whether or not the (/)-theory 
Eq. 12 exhibits off-diagonal long-range order and a bro- 
ken t7(l)-symmetry. An entirely equivalent interpreta- 
tion of the change in Ol was given in Ref. which did 
not involve a local field 0(r), but number conservation of 
vortex-lines threading the entire superconductor. This 
number is conjugate to the phase-field of the local com- 
plex field (j){r) . An advantage of the present formulation 
involving Eq. |l^ is that it directly relates he change in Ol 
to the long-distance part of a correlation for a local field, 
and hence to a local order parameter < (j>(r) >. This 
connection makes it at least plausible that the change in 
vortex-tangle connectivity, i.e. a change in the geometry 
of the vortex-tangle, may be related to a thermodynamic 
phase-transition. We emphasize that the present problem 
is very different from the percolation transition known to 
occur in the 3D Ising-model, and which has nothing toJxi 
with the thermodynamic phase-transition in that modeEB, 



8. Structure function S(k) 

To probe the structure of the VLL, we consider the 
structure funC|tijm,for vortex segments directed along the 
field directionLJ'Ej. Given an applied field along z-axis, 
the structure function 5(k) is defined by 



5(k) 



<l Sr"-^('')6xp [»k-r] p> 
{fVf 



(24) 



Here, k is a reciprocal lattice vector 
fraction / is defined in Eq. 



26 



and the filling 
In the VLL phase, 
5(kj_,fc2 — 0) has (5- function Bragg peaks at k_L = K, 
where K is reciprocal lattice vector for the ordered VLL. 
When the VLL melts, the Bragg peaks are smeared out 
5(K, kz = Q) drops to a very small value. Thus, the VLL 
melting temperature can be defined as the temperature 
where S'(K, kz = 0) shows a sharp drop to a value very 
close to zero. In a liquid where we have full rotational 
invariance, we expect 5(k^, kz = 0) to exhibit a ring pat- 
tern, with a very small amplitude. 
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9. Extended Landau gauge: 

Periodic boundary conditions together with Landau 
gauge 

Ay = 2TTfx, 

give rise to a constraint, L^f ~ 1,2,3... Thus, for given 
Lx , the smallest / allowed is / = 1/L:r . To perform sim- 
ulations and finite size scaling of systems with very low 
filling fractions, we define an "extended" Landau gauge, 

A = ^-f ; Ay^ , (25) 

J^xJ^y i^xJ^y 

where Ux, n, my, m are positive integers satisfying Ux n = 
Ly, and viy m = Lx- The filling fraction / is now given 
by 

^ _ n m [ux - my] 

LxLy 

Hence, it is possible to chopse systems with a filling frac- 
tion as low as / = l/LxLj^. 



B. Details of the Monte-Carlo simulations 

The statistical mechanics of the 3DXY model and the 
FG model is investigated by Monte-Carlo simulations on 
the effective Hamiltonians Eq. |l^ and Eq. 0. 

For the 3DXY model, a Monte-Carlo move is an at- 
tempt to replace a phase angle at a given site 9{r) with a 
new randomly chosen phase angle 0' S [0, 27r >. For the 
FG model, a Monte-Carlo move is an attempt to replace a 
complex number at a given site ijji^) with a new randomly 
chosen complex number V''- Here, \ip'\ G [0 : 1 -I- e] and 
6' G [0, 27r >. We have introduced a small positive pa- 
rameter e to allow the system to perform Gaussian fluctu- 
ations, around the extremal field configuration — 1, 
at very low temperature. Note that we are letting the 
amplitude fluctuate around its mean value at every tem- 
perature. The Monte-Carlo move is accepted owejected 
according to the standard Metropolis algorithmEJ. If the 
new phase angle causes a gauge invariant phase differ- 
ences jfj_ = Af^0 — Afj, to fall outside the primary inter- 
val [— TT, TT >, we take it back into the primary interval. 
This compactization procedure creates a closed unit vor- 
tex loop around the link where is changed. In this 
way, all the vortex loops introduced into the system are 
closed, and the net induction is always conserved. 

A Monte-Carlo sweep consists of Lx x Ly x Lz Monte- 
Carlo moves. Typical runs consist of 1.2 x 10^ sweeps 
per temperature, where the first 2 x 10"* sweeps are dis- 
carded for equilibration. Near the phase transitions up 
to 2 X 10^ sweeps per temperature is necessary to capture 
the correct physics. For a given system, we always start 
the simulation by a cooling sequence, where the starting 



temperature is significantly higher than all temperatures 
associated with phase-transitions or crossovers the model 
might exhibit. The results shown in this paper originate 
both from cooling and heating sequences. Since these 
two methods give essential identical results, we do not 
differentiate between them. 

In order to resolve anomalies in the specific heat, we 
must in some cases perform simulations on systems as 
large as 360'^. To be able to carry out simulation on such 
large systems, we must 1) write part of the code in assem- 
bly and 2) carry out the simulations in a parallel manner. 
Our systems are divided into "black and white" subsys- 
tems, arranged in a 3D checkerboard pattern. Each black 
subsystem has only six white subsystems as its nearest 
neighbors, and visa versa. Since the 3DXY and the FG 
model only have nearest neighbor interactions, all subsys- 
tems with the same color can be updated simultaneously. 

To be able to calculate a nonparallel-able routine as Ol 
in an effective manner, we divide the computer nodes in 
2 groups, the large main group takes care of the Monte 
Carlo simulation, and a small subgroup carries out, si- 
multaneously, the calculation of O^. 

C. Model parameters 

a. System sizes: We put our coordinate (x,y,z)-axes 
along the crystal (a,b,c)-axes. For the anisotropic cases, 
we assume uniaxial anisotropy, and use the crystal c- 
axis as the anisotropy axis. We perform simulations 
on tetragonal systems with dimensions Lx, Ly, Lz- The 
main part of the simulations is done on cubic or nearly 
cubic systems. Nearly cubic systems Lx ^ Ly — Lz is 
some times necessary in order to satisfy the boundary 
conditions enforce by the extended Landau gauge, Eq. 
p5| . To check for the finite size effect of Ol, we carry 
out simulations on slab systems with the aspect ratios 
Ly/Lz = 1.00,1.25,1.5,1.75,2.00. System sizes up to 
360'^ were used. 

b. Cooper-pair chemical potential a{T): . We let 
the Cooper-pair chemical potential have the simple linear 
form 

a{T) _ T - Tmf 
a{0) Tmf 

We have also carried out simulations with other forms 
for a(r)/a(0), such as tanh([T - Tmf/Tmf]To). Here, 
To is a constant regulating the size of the region where 
a{T)/a{0) grows from -1 to 1. The results are, how- 
ever, qualitatively the same as for the linear case. The 
parameter Tmf is the parameter effectively controlling 
the width of the critical region in these calculations. In 
units of Jo, Eq. ||, we write Tj^jp = kBTMF / Jo- The 
values we will use are T'^jp = 0.3,1.0. An estimate for 
what temperatures these values correpond to may read- 
ily be obtained by using = 11 A, \ab = 1500 A, im- 
plying that T'f^,jp = 1 corresponds to 300isr, while us- 
ing \ab — 2000 A implies that T^jp = 1 corresponds to 
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180-fs'. These are very reasonable numbers. With these 
parameters, it will be shown that amplitude fluctuations 
of the local Ginzburg- Landau order parameter, when in- 
cluded on an equal footing with the phase-fluctuations, 
are far from being critical. 

c. Anisotropy parameter T : The anisotropy parame- 
ter r is defined as 

Note that F > 1 only when the layering of the supercon- 
ductor to be simulated is pronounced, i.e. > for at 
least one direction fi. In this article, we consider systems 
with the anisotropy parameter F = 1, 3, 7. 

d. Filling fraction f : The filling fraction along the 
/i-direction, is defined as 

2nf^ = (A X ^)^. 

is a measure of the fraction of fiux quanta of magnetic 
induction penetrating a single plaquette with surface nor- 
mal along //. When the magnetic field is applied along 
z-axis, f^^fy^ 0, and 

Ba^ay ,DXY BC, . 

In this work, we consider filling fractions / = 
0,1/20, ,1/1560. 

IV. MONTE-CARLO SIMULATIONS, B = 

In this subsection, we discuss the zero field 
superconducting-normal (S-N) phase transition, both in 
terms of the usual Ginzburg-Landau order field V'(r), and 
in terms of the behavior of topological excitations which 
can be tied to the formulation of the transition using the 
disorder- field picture presented in Section IIF. We com- 
pare our results obtained from the FCj model to known 
simulation results of the 3DX"*>pjpodeOE3, the London 
modelEl, and the Villain modelElo. 

Unless otherwise stated, in this subsection we show 
simulation results for the FG model with the parame- 
ters: / = 0, F = 1, = 0.3,1, a^/e,, = 6, and 
V = 60^. We have chosen a^/^^ — 6 to slightly enhance 
the critical features of the FG model. Simulations of the 
FG model using a smaller ratio a^/^^ ~ 4 leads to the 
same conclusions, but larger systems and longer simu- 
lation times are required to obtain the same quality of 
the data. What we will find is that the width of the re- 
gions where phase-fluctuations dominate is controlled by 
the parameter Tmf, increasing with Tmf- Our picture 
of the zero-field transition as a vortex-loop unbinding is 
however borne out regardless of the value of Tmf, close 
enough to the true critical temperature T^. 



A. Order field 

First we present results for the S-N phase transition 
in terms of an order field picture, i.e. in terms of the 
ordinary local pair- wave function ip{r) of the Ginzburg- 
Landau theory. We obtain results of the full GL the- 
ory including amplitude and phase-fluctuations. What 
will be shown, even when amplitude fluctuations are in- 
cluded, is that we obtain a clear picture of the S-N phase 
transitionJiiilerms of an Onsager-Feynman vortex loop 
unbindingfju, driven exclusively by topological phase- 
fluctuations of the Ginzburg-Landau order parameter. 

In Fig. we plot the helicity modulus Y^, the local 
density of Copper pairs < 1-0' P >, and the superfluid 
condensate density | <■(/''> P as functions of tempera- 
ture. We see that the condensate density | < -0' > p is 
zero above a critical temperature and develops a finite ex- 
pectation value below Tc- In contrast to this, < l^'p > is 
finite both above and below Tc- Close to T = Tc we have 
performed the simulations for a very dense set of temper- 
atures, and from the top panel of Fig. 0, we may however 
discern a kink in the curve and hence a singular behavior 
of the temperature derivative of < |0'p >. The top panel 
shows the results for Tmf — 1-0, while the lower panel 
shows the same results for Tmf = 0.3. The difference be- 
tween the two panels is that since Tmf has been changed, 
the width of the critical region has changed, increasing 
upon increasing Tmf- Had we chosen Tmf = 0.01, an 
appropriate value for conventional superconductors, the 
curves for | < ip' > p and < l^'p > would have been 
indistinguishable, the conventional BCS mean-field pic- 
ture of the superconducting transition would have been 
appropriate. The reason that it is no longer the case in 
the high-Tc cuprates is the large energy scale for pairing, 
coupled with the fact that the phase-stiffness is low. 

Note how the curve for ~ |t-|2/3-i)i' bends slightly 
more sharply towards zero than the curve for the conden- 
sate density \ < ip' > l"^ ^ hl'^'^ : as expected for a positive 
rj, since in that case 2(3 — rjv < 2(3. In fact, this provides 
a nice consistency-check on the Monte-Carlo simulations. 

We may therefore conclude that the existence of 
Cooper pairs does not imply superconductivity, and the 
local Cooper pair density can not be used as a probe of 
the superconducting phase. Rather, the superconducting 
phase is characterized by a finite expectation value of the 
superfluid condensate density, i.e < 0'(r) >^ 0. Recall 
that I < 0' > |2 ^ only if < V'(r) >^ 0. The slow 
decay of < |V''P> above the true superconducting tran- 
sition temperature Tc is due to gaussian fluctuations in 
the amplitude of the order parameter. 

The results shown in Fig. |l| justify, quite clearly, that 
the neglect of amplitude fluctuations in order to study 
the critical behavior, is an entirely appropriate approach 
to this problem. In the normal phase, the phase angle 
of the order field is uniformly distributed. Fig. |l| inset, 
while for T < T^, the system spontaneously chooses a 
preferred phase angle giving a peak in De{9). Due to 
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our finite set of discrete phase angles, Dg{6) = l/Ng 
T > Tr. and not zero as in the continuum-0 limit. 



for 




FIG. ^ Helicity modulus Tz, local Cooper-pair density 
< It/j'I^ >, and superfluid condensate density | < > p as 
functions of temperature for the Ginzburg-Landau model in 
a frozen gauge approximation. Upper panel shows results for 
/ = 0, r = 1, T'mp = 1.0, Op/^^ = 6, and V = 60^ Lines are 
guide to the eye. Tz and | <ip' > P develop finite expectation 
values for T < Tc, while < l?/)'!^ > is finite both above and 
below Tc- Inset: The distribution function of the phase angle 
of the order field Dg{9) as a function of 9, for several tem- 
peratures. Below Tc, a preferred phase angle is chosen and 
the global U{1) symmetry is spontaneously broken. Note how 
the phase fluctuates around a mean value even in the ordered 
phase. Lower panel shows the same for T'l^ip = 0.3. The 
width of the critical region has decreased, but the amplitude 
fluctuations are still far from being critical. 



Note the difference between this picture compared to 
the usual mean field picture of the S-N phase transition. 
In the mean field picture, the phase angle of the order 
field does not fluctuate. Thus, < IV^'P >= \ < i/j' > 



and Tc can be defined as the highest temperature where 
the local Cooper pair density < j-^'P > develops a fi- 
nite expectation value. This happens, at the mean field 
level neglecting all fiuctuations in ?/)', at T]\if where the 
Ginzburg-Landau parameter a{T) = 1 — T' /T'j^jp be- 
comes negative. Including fluctuations < >y^ \ < 
ip' > and Tc is defined as the lowest temperature 
where the superfluid condensate density | <-!/;'> ^ still 
maintains a value of zero. 

Given that the condensate density is non-zero be- 
low Tc, we next focus on a global quantity, the long- 
wavelength limit of the helicity modulus T^, or equiva- 
lently the superfluid stiffness in the /^-direction. In Fig. 
|l| we see that vanishes for temperatures T >Tc, and 
develops an expectation value for T < Tc- Thus, the 
superconducting phase exhibits global phase coherence, 
while the normal phase does not. We have also calcu- 
lated Ta; and Tj,, and found (not shown) that they show 
the same behavior as T^. Apart from minor details, we 
see in Fig. ^ that the helicity modulus is proportional 
to the condensate densityO. We will also show that this 
equality also applies to the finite fleld case. 

At low temperature, T^^ decreases linearly. This fea- 
ture is also obtained in the zero field— 3DXY modeH, 
but not in the zero field Villain modeltZl. In the Villain 
model, the spin waves and the vortex loops can be ana- 
lytically decoupled. Here, at low temperatures, spin wave 
excitations do not affect the vortex loops excitations and 
the superfluid phase stiffness should decay in an acti- 
vated manner due to the excitation of vortex loops. In 
the 3DXY model, the spin wave and the vortex loops are 
coupled together. Whether or not the low-temperature 
features of in Fig. ^ can explain experimental data 
on the temperature dependence of 1/A^, see for instance 
Rcf. ^ is an interestinff-jlput-ao far unsettled issue, see 
also the results of RefsJl3HE3. Within the anisotropic 
3Z)Ari"-model, the helicity modulus has a larger, but 
negative slope of its linear low-T behavior compared to 
Ta; and Ty. On the other hand, it is not entirely trivial 
to connect T z{Tl)JT to the T — 0-normalized super- 
fluid density Pszl3, which iS|thc quantity measured in the 
experiments of Hardy et alsB. However, our main point 
of emphasis is that the vanishing of the superconduct- 
ing phase stiffness at T — Tc is caused exclusively by an 
unbinding of large vortex loops. Further evidence for the 
connection between and the vortex loops can be found 
in simulations of the lattice London model- ssihere vortex- 
loops are the only degrees of freedon£illi3'E3. Here, the 
normalized helicity modulus limfc^o Tf^(T)/T^(r = 0) is 
renormalized to zero at Tc exclusively by the expansion 
of vortex loops. 

The critical behavior of the 3DAry-model including 
coupling between spin-waves and vortices, is the same 
as for the vortex-content of the same theory, but taken 
in the Villain-approximation, where no coupling between 
vortices and spin- waves exist. The only excitations in the 
C/(l)-symmetric gauge-theories in the frozen gauge- and 
amlpitude-approximation, are longitudinal and trans- 
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verse phase-fluctuations, i.e. spin-waves and vortex- 
loops. 

The above considerations and results provide an over- 
whelming amount of evidence in favor of the proposition 
that unbound vortex-loops are precisely the critical fluc- 
tuations of an extreme type-II superconductor. 

Approaching Tc from below, T^(T) decays to zero 
with an exponent consistent with v — 2P — rju, and 
77 « 0.04 as discussed in Section III A. 5, see Fig. ||. 
For the special case of d = 3, which we consider, we have 
Ps/i ~ ~ kr, and hence we find 2(3 — rjv = v < 2(3. 
We will show in the next section that this renormalization 
of is nicely explained by the expansion and blowout 
of thermally excited vortex loops. 



bottom panel shows the vortex-loop distribution func- 
tion D{p) as a function of perimeter p for a number of 
temperatures T < Tc, while the inset of the bottom 
panel shows the temperature dependence of the long- 
wavelength vortex-line tension s(T). In the top panel 
it is clear that the loss of helicity modulus, the anomaly 
in specific heat, and the abrupt change in Ol all occur 
at precisely the same temperature. The change in the 
decay of D{p) also occurs at the same temperature, Tc. 



B. Topological excitations and vortex-line tension 



We now discuss the S-N phase transition in terms of 
the behavior of topological objects of the FG model, i.e. 
closed vortex loops and vortex lines. 

The first figure of this section shows our probe of the 
vortex-tangle connectivity, Ol, as a function of tempera- 
ture, for various system sizes and F = 1 and 7. We reem- 
phasize that Ol is not suggested as an order parameter 
of the transition, but that it nevertheless can be tied to 
a local order parameter via the formulation of the theory 
given in Section IIF, unlike the [known spin-percolation 
transition in the 31? Ising-modeO. 

In Fig. |[ we show the probability of finding a con- 
nected vortex-tangle across the system in zero magnetic 
field, for F = 1, 7, and system sizes L^, with L — 6, ..64. 
Notice how the curves cross at approximately the same 
temperature and get progressively sharper. Similar re- 
sults were seen for considerably smaller system sizes 
L = 4, 6, 8 in Ref. Below, we will also give re- 

sults for much larger system sizes, confirming that the 
crossing temperature in Fig. || gives a good estimate 
for the threshold temperature for vortex-loop unbinding 
throughout the sample. As pointed out in Ref. |9|, such a 
finite-size effect indicates that a percolation threshold 
ists for the vortex tangle in the thermodynamic limitE3. 
A precisely similar finite-size effect in Ol will be seen 
in finite magnetic field, to be considered in Section VB. 
This will happen inside the vortex liquid phase at ele- 
vated magnetic fields, but will coincide with VLL melting 
at low fields, and suggests the revision of the picture of 
the molten phase of the Abrikosov vortex system purely 
in terms of a vortex-/me liquid. 

We next proceed to correlate the change in Ol with the 
unbinding of large- vortex loops and the loss of vortex- line 
tension at T = Tc, by correlating its abrupt change with 
the characteristics developing in D{p), which probes the 
typical size of thermally induced vortex loops in the sys- 
tem. We first consider the case of the 3£'Xy-model, for 
which the results are shown in Fig. |^. The top panel 
shows specific heat, Ol and helicity modulus, while the 
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FIG. ^ Ol as a function of temperature for several sys- 
tem sizes for the 3DXY model with f=0. Lines are guide to 
the eye. Top panel: F = 1, bottom panel: F = 7. Note 
the finite-size effect in Ol, with the crossings of the curve 
approximately at the same temperature. 
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of the zero temperature fixed point. At and above Tc, 
Ol = 1, and vortex loops with infinite size always exist. 
The length scale Lq{T) has diverged, showing that there 
are vortex-loops on all length-scales with a power-law tail 
in the distribution. Such loops cannot be coarse grained 
away and taken into account by any "appropriate renor- 
malization" of the zero-temperature theory. Thus, the 
S-N phase transition can be viewed as a blowout out of 
thermally induced vortex loops. Above Tc, free thermally 
induced "vortex lines" exist in all direction, and any in- 
finitesimal applied current will move thesCr-thermally in- 
duced "vortex lines" and dissipate energytll. Thus, the 
system is in the normal phase. 
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FIG. ^ Top panel: Specific lieat C, helicity modulus Tz, 
and Ol as functions of temperature for tlie 3DXY model with 
/ = 0, r = 1 and V = 120'^. Lines are guide to the eye. 
Bottom panel: Vortex- loop distribution function D{p) as a 
function of loop-perimeter p for various temperatures. Lines 
are fits using D{p) = p"^/^ exp{-~£{T)p/kBT). At T = Tc the 
decay changes from exponential to algebraic implying that the 
vortex-line tension e vanishes. Inset of bottom panel shows 
e{T). Solid line is a fit using |T' - T^'^, with 7 = 1.45 ± 0.05. 



For T < Tc, Ol = 0, and all vortex loops are confined, 
with typical size given by i^oiT) = ksT / e{T), where 
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where s{T) is the effective long- wavelength vortex line 
tension, equivalently the free energy per unit length of 
vortex-lines. These objects, present also in the low- 
temperature phase, cause only a local perturbation of 
the order parameter in the system, and may simply be 
"coarse grained" away. The low energy physics of the 
model is therefore described essentially by the physics 



1e-06 



1e-07 



1e-08 



FIG. Q. Top panel: Specific heat C, helicity moduli 
T^jTz, and Ol as functions of temperature for the 3DXY 
model with / = 0, F = 7, and V = 140^. Lines are guide to 
the eye. 

Bottom panel: Vortex-loop distribution function D{p) as a 
function of loop-perimeter p at various temperatures. Lines 
are fits using D{p) = 0.37p"^-^^ exp(-e(r)p/fcsr). Inset: 
Vortex-line tension e(T) as a function of temperature. Solid 
line is a fit using \T' - T^p , with 7 = 1.45 ± 0.05. 
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the sharpness of the manner in which the moduh ap- 
proach zero at Tc, there is no high-temperature tail as 
one would have found in too small systems. This in fact 
serves as a highly non-trivial benchmark on the quality 
of the Monte-Carlo simulations. 
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FIG. 1^. Top panel: Specific heat C calculated using 
Eq. |l^, Ol, helicity modulus T^, and superfluid density 
I < Tp' > P as functions of temperature for the Ginzburg- 
Landau model in a frozen gauge approximation with f — 0, 



1, T' 



1, a^/^M 



6, and size V — 60 . Lines are 



guide to the eye. Inset shows the specific heat calculated us- 
ing Eq. 

Bottom panel: Distribution of vortex loops D{p) as a function 
of the loop-perimeter p for several temperatures. Lines are fits 
using = L15p~^/^ exp(-e(T)p/fcsr). Inset: Vortex-line 
tension e{T) as a function of temperature. Dotted line is a fit 
using \T' - T^p, with 7 = 1.45 ± 0.05. 

In Fig. ^, we show the specific heat anomaly, the he- 
licity moduli Tx and T^, as well as Ol for the 3DXY- 
model, with F = 7. The correlation noted above in con- 
nection with Fig. ^ is again perfect, the only difference 
being that the specific heat anomaly has become more 
symmetric due to the increased anisotropy, F = 7. Al- 
though the amplitude of T^, is larger than the amplitude 
of Tz due to the uniaxial anisotropy along the z-axis, 
the temperature at which they vanish, and the power 
law with which they vanish, are the same. Note also 
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Figs. H and ^ show essentially the same as Fig. ^, but 
now for the FG model, i.e. including amplitude fluctu- 
ations on an equal footing with the phase-fluctuations. 
Clearly, the picture that it is the topological phase- 
fluctuations, or the vortex-loop unbinding, that drives 
the superconductor normal-fluid transition, is not at all 
altered by the fact that amplitude fluctuations are in- 
cluded. This is a reconfirmation of the results obtained 
in Section IVA, showing that amplitude fluctuations of 
the local Ginzburg-Landau order parameter have a large 
mass at the critical temperature where the superfluid 
density vanishes. The vortex-loop unbinding, which is the 
microscopic mechanism driving the transition, cannot be 
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reparametrized in terms of critical amplitude fluctuations 
of the local order parameter of the theory, since the latter 
are clearly nowhere close to being critical. 

For a more detailed study of the properties of thermally 
induced vortex loops, we now focus on the vortex-loop 
distribution function D{p) as a function of vortex-loop 
perimeter p at various temperatures. These are shown 
in Figs. 1^, 1^, and ^ and are clearly well approximated 
by the form D{p) = Ap^^/^ ^-<^{T)p/kBT j^j. g^jj tempera- 
tures considered. Note that s(T) is the only temperature- 
dependent fitting parameter in all plots. The effective 
long-wavelength linetension of vortex loops is finite be- 
low Tc, and vanishes for T > Tc. The physical picture 
of this phase transition is as follows. Below Tc, e{T) is 
finite defining a typical length scale for the vortex loops, 
Lq = ksT / e{T). Here, Diji) is dominated by an expo- 
nential decay and vortex loops with much larger perime- 
ter p than Lqj ^re exponentially suppressed. Thus, the 
topological excitations that are present in the system may 
be coarse grained away. At and above Tc, e{T) = and 
no typical length scale for the vortex loops exist; the 
length scale Lq has diverged. Here, D{p) is purely al- 
gebraic, and vortex loops of all sizes including infinite 
size, exist. Thus, the S-N phase transition at Tc is trig- 
gered by an unbinding of large vorte* loops, analagous 
to the Onsager-Feynman mechanism,tZI suggested for the 
superfluid-normalfluid transition in *_ffeE3. 

In the insets in the bottom panels of Figs. |^, ^, |^, and 
^, we show the vortex line tension e{T) extracted from 
the vortex loop distribution function D{p). Regardless of 
whether the iDXY- or the FG-models are used, we find 
that the long-wavelength vortex-line tension vanishes as 

e(T)-|T-Tcp; 7 = 1.45 ±0.05. (28) 

The numerical value of the exponent 7 has been extracted 
from the systems with the largest critical regions, i.e. 
Figs. ^, 1^, and ||. The system shown in Fig. |^ does not 
allow a very precise value for 7 to be obtained, although 
the qualitative aspects of the results are clearly precisely 
the same as those for the 'iDXY -model and the FG- 
approximation of the GL- model with Tjjj^ — 1.0. This 
implies that the typical vortex-loop perimeter diverges 
when Tc is approached from below, using Eq. as 

To(T)~ |T-Tc|-^ (29) 

such that Lq{T) is a power of the correlation length ^ of 
the 3DXY model. 



C. Anomalous dimension of the dual field 

We next connect the result for e{T) to the anomalous 
dimension of the dual field (j). It is natural, within the 
formulation of the problem given in Section IIF, to as- 
sociate the proliferation of unbound vortex-loops with 



a vortex-loop susceptibility, or equivalently a suscepti- 
bility for the 0-field of Section IIF. This is seen as fol- 
lows. The proliferation of unbounded vortex loops as the 
temperature of the superconductor is increased, is asso- 
ciated with the development of long-range correlations 
in the two-point correlation function of the dual field, 
G{x) =< 0*(a;)(/>(O) >, where on the low-temperature 
side the dual order parameter has zero expectation value, 
< (j) >— 0. A scaling Ansatz for G{x) reads 

G(^) - ^T^T^ Si-/0, (30) 

where rj^ is the anomalous dimension of the dual field 
4>, ^ is its correlation length, and G{x/£,) is some scaling 
function. The square of the mass of the dual field, m^, 
is therefore naturally mapped to the line-tension e(T) of 
the vortex-loops. This follows from the observation that 
the dual boson system of which the (/)-theory is a field- 
theory description, has a chemical potential which in 
turn is nothing but the line-tension s{T) of the vortex- 
loop system, when the density distribution D{p) is viewed 
as a partial density in a fugapty expansion for the den- 
sity of the dual bose-systemO. The Fourier-transform 
G{k) (p* {k)(j>{~k) > of G{x) may be written on form 

G{k)^e-"'' HK), (31) 

where !F{ki) is some new scaling function. The k —t Q 
limit of this is the static uniform susceptibility X4> of the 
dual field on the low-temperature side, where < (p >= 0- 
On the other hand, as long as the dual field is massive, 
which it is on the low-temperature side, we must have 
limfc^o G{k) = rn'^'^ . Hence, we obtain 

XA,r^\^-^ e^-"* - -"^42-^,^). (32) 
e 

The field (f> has a correlation length exponent rgiven by 
— 2/3, the same as for the SUATF-modelO This 
follows from the fact that it is a thermodynamic ex- 
ponent describing the divergence of one and the same 
length in the Ginzburg-Landau theory and dual theory. 
Very importantly, it must be equal both for the dual 
model aad-its Ginzburg-Landau counterpart by "strong" 
dualityE3~EJ. Were this not to hold, the dual of the dual 
theory would not be the original theory, as it ought to 
be. The aboye of course precisely amounts to the Fisher 
scaling-lawEll relating the susceptibility exponent of the 
dual field 7^ to and r\^ 

70 = (2- J70). (33) 

Using our estimate 7^ — 1.45 ± 0.05 with = 2/3 gives 
T]^ = — 0.18=f0.07 in close agreement with previous renor- 
malization group calculationso, who found r\^ = —0.20 
to one-loop order. 

The result r\^ — —0.18 =F 0.07 obtained directly from 
computing the statistics of the loop-excitations of the 
31) ATF- model is a truly noteworthy result, when viewed 
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juxtaposed to the RG-calculations of Ref. In Ref. |9^, 
the RG-result for the anomalous dimension of the dual 
field was obtained directly from the dual theory. On the 
other hand, our numerical result is obtained directly from 
the phase-only approximation to the original Ginzburg- 
Landau theory. The agreement shows conclusively, and 
to our knowledge for the first time, that viewing the zero- 
field transition of the 3D Ginzburg-Landau theory as a 
vortex-loop unbinding, which is the phase-transition of 
the dual theory, is pre(ds.ely correct, not only qualita- 
tively, but quantitativel'^3. 

At and below Tc, the order field < ip'i^) > develops an 
expectation value, and explicitly breaks the global C/(l) 
symmetry of the GL-theory. In contrast to the order 
field picture, in a description using only topological exci- 
tations, the global C/(l) symmetry is hidden. There does 
not appear to be any symmetry operation involving the 
phase of a local field, that will leave the effective action 
Eq. |l^ invariant. Therefore, there is also no obvious 
local quantity that develops an expectation value in the 
non-symmetric phase. Nevertheless, it is possible to de- 
fine a global quantity that implicitly probes the breaking 
of the global f/(l) symmetry, namely Ol- Let denote 
the number of "vortex lines" (percolating directed vortex 
paths without using PBC) along the /i-direction. Below 
Tc, Nf^ is fixed to zero and Ol = 0. Concomitant with 
the conservation of the global quantity N^^, the system 
must exhibit a global U{1) symmetry. At and above Tc, 
iV^ develops an expectation value and Ol 7^ 0. This 
leads to a broken C/(l) symmetry. 



V. MONTE-CARLO SIMULATIONS, B / 

We next discuss the indications we have of phase- 
transitions in the vortex-system in a finite magnetic field. 
In addition to the first order VLL melting transition line 
T„i{B) which we map out for a large range of filling frac- 
tions, we find indications for a new phase transition in the 
vortex liquid. We emphasize that in all simulations per- 
formed in finite magnetic field, the filling fraction is low 
enough to ensure that there is zero transverse Meissner 
effect at any temperature of interest. That is to say, the 
vortex-line lattice is depinned from the numerical lattice 
at much lower temperature than the temperatures where 
the Bragg-peaks in the structure function of the VLL van- 
ishes. Therefore, commensuration effects due to defining 
the theory on a lattice effectively have been eliminated 
at the temperatures of interest. 

Before entering into the discussion, a clarifying remark 
is appropriate. Note that the phase-transition that we 
suggest may he taking place inside the vortex liquid, is 
not a transition from a disentangled low-temperature vor- 
tex liquid to an entangled high-temperature vortex liquid, 
as discussed by numerous previous authors. Such a tran- 
sition would have has its hallmark that the superfluid 
stiffness along the magnetic field, or equivalently the he- 



licity modulus T^, would vanish inside the vortex- liquid 
phase. This ha«,iioj«-been conclusively demonstrated not 
to be the casetaiilN. 

At the phase transition we propose inside the vortex 
liquid, a global U{1) symmetry is spontaneously broken. 
The probe for this new phase transition is, as for the zero 
field case, a sudden change in the connectivity of the vor- 
tex system, probed by Ol- As defined in the zero-field 
section, we have found that the vortex-loop distribution 
function D{p) is not useful as a probe for the vortex line 
tension in finite field. Fortunately, non-local quantities 
with directions built-in as the helicity modulus, and the 
percolation probability still apply in finite field. To probe 
the VLL melting, we must define a new direction depen- 
dent non-local quantity, namely the structure function 
for vortex segments parallel to the field, S'(k). Quanti- 
ties probing local properties as < > and | < V'' > | 
are still useful in finite field. 



A. First order VLL melting transition 



In this sub-subsection, we show and discuss our results 
with the focus on the first-order VLL melting transition. 

We show in Fig. |^ the structure function for vortex seg- 
ments along the field direction 5'(K = [tt/IS, 47r/15, 0]) 
and the fluctuation specific heat C as functions of tem- 
perature for the FG model. Here, K is a reciprocal lattice 
vector of the VLL. 

We see in Fig. ^ that S'(K) shows a sharp drop at 
ksTm/Jo = 0.26 from ~ 0.2 to 0, indicating a flrst or- 
der VLL melting transition. The Lindemann number for 
this melting transition estimated from the Debye- Waller 
factor at the melting temperature, is 0.25, consistent 
with previous estimatesB^EZHlj. This abrupt disappear- 
ance of the VLL structure jfunction has also been experi- 
mentally found in BSCCOl23a. In the high-field regime, 
the position of the ru^liing line is well estimated by the 
Lindemann- criterioru'Cn. 

Density plots (not shown) of S'(k_L,0) show Bragg 
spots for T below T^, and rings pattern for T just above 
Tm. Here, = [k^^ky^O] and (k^^ky) G [-tt : tt]. 
This expected peak-to-ring featurCpj* also found at Tm 
m sirnulations on the 3DXY modejlj and in the Villain 
modelEj'EZl. Below T,„, the vortex lines perform small 
Gaussian fluctuations around their equilibrium positions, 
and rcnormalizc 5(K; T) to a smaller value compared to 
S(K;T = 0). Precisely at the speciflc heat shows a 
delta-function like anomaly. Fig. ^ indicating the well- 
established flrst-order character of the VLL-melting tran- 
sition in clean systemsEilEl. Note that T,„/Tbc2 increases 
for decreasing Tmf- For very small Tmf, where the mean 
field approximation applies, Tm — Tbc2- 
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FIG. 1^ Top panel: Structure function for vortex segments 
along the field direction ^(K = [tt/IS, 47r/15, 0]) and specific 
heat C as functions of temperature for the Ginzburg-Landau 
model in a frozen gauge approximation with / = 1/60, F = 7, 
T'mf ~ Ij 'Jm/Cm = 6, and V = 60^. ^(K) shows a sharp 
drop at kBTm/ Jo =0.26 indicating a first order VLL melting 
transition. C shows a delta-function like anomaly at Tm, the 
VLL melting transition, detailed in the inset. 
Bottom panel: Same as top panel, but for T^fp = 0.3. 

In Fig. ^, we show the local Cooper-pair density 
< 1-0' P >, the local condensate density | < tp' > p and 
the hclicity modulus along the field direction as func- 
tions of temperature for the FG model. 

We see in Fig. ^ that < \ip'\'^ > does not show any 
particular anomalous feature at the melting tempera- 
ture Trn, it is completely smooth. In contrast to this, 
I <ip' > p shows a sharp drop to zero at T™, which in- 
dicates a first order phase transition at T^- < I'^'P > 
shows a small tail above T^. However, note that this tail 
decreases towards zero for increasing simulation length. 
Thus, the superfluid condensate density is finite in the 
VLL phase, jumps dis continuously to zero at Tm, and is 



zero in the entire vortex liquid phase. The vortex-liquid 
phase is therefore phase-incoherent, and we denote the 
phase as an incoherent vortex liquid. The behavior of the 
helicity modulus shows that there is no signature of any 
transition from a disentangled to an entangled .vnrtex liq- 
uid phase, in agreement with previous resultsHu. In the 
language of the non-relativistic 2D boson-analogy com- 
monly invoked in the study of vortex-liquids, there is no 
transition from a normal fluid to superfiuid bose-system 
at T = of the bose-system, in agreement with general 
arguments given by Landau. In the absence of disorder, 
a normal bosc-fluid cannot exist at T = in the thermo- 
dynamic limit, the only possible phases are bose-crystals 
or superfluids and hence there should be a direct tran- 
sition from an insulating bose-system (corresponding to 
the Abrikosov vortex lattice phase) to a superfluid bose- 
system, corresponding to an incoherent vortex liquid. 

This is precisely what is observed in all our simula- 
tions, and demonstrates that all our systems are suffi- 
ciently large in the z-direction to capture this physics cor- 
rectly. In Fig. H the helicity modulus along the field di- 
rection T z also shows a sharp drop to zero at Tm. Thus, 
the vortex liquid is incoherent and can not carry any su- 
percurrent. Note that is a global quantity and says 
nothing about local (in space and time) phase coherence. 
Thus, even when = 0, the system can exhibit local 
superconductivity, and diamagnetic response is expected 
in the incoherent vortex liquid. 

It is interesting to note that while this comes out 
fairly |-s|p|ightforwardly in simulations on the 2>DXY 
modelllj'E^, it is much more difficult to obtain us- 
ing Quantum Monte Carlo simulations of the 2D non- 
relativistic boson-analogv-jjif, the lines-only approxima- 
tion to the vortex-liquidoEj. This is due to the fact 
that the so-called winding number, which is the appro- 
priate quantity to measure, is essentially inaccessible for 
large systems. One has to resort to computing quantities 
which are not the same as the winding number, but hop- 
ing that it is a representative of it. For a nice discussion 
of this point, see Refs. It still has not been very 

well established precisely which fluctuations are respon- 
sible for destroying longitudinal phase-coherence at the 
VLL melting transition. 

The loss of global phase coherence in the incoherent 
vortex liquid does not mean that the layers are decou- 
pled. The correlation length ^^(T), for all directions /i, 
jumps from infinity to a finite value at Tm, and further 
decreases for increasing temperature. Estimates for var- 
ious correlation lengths along the field direction can be 
found in Refs. [i^JOq . This picture nicely explfiins the ex- 
perimental results found in BSCCO in Ref.O. Here, the 
difference in the high frequency resistivity, for different 
sample configurations, is used to probe the VLL melting 
transition. In the VLL phase is infinite in all direc- 
tions, and the "rigid" VLL does not contribute to the 
dissipation when a higltfrequency circulating current is 
applied (Fig. 2 in Ref.E3) is applied. In the vortex liq- 
uid however, this dissipation depends on the correlation 
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length in the different directions. Due to anisotropy and 
the appUed field along the crystal's c-axis, 7^ (,ab, and 
the high-frequency microwave dissipation is different for 
different sample orientation. For increasing temperature, 
^^(T) decreases further ^-and thus the dissipation rate in- 
creases as found in Ref.O. 



ence along the field direction is still intact, and does not 
disappear before a new "entangled" phase transition at 
Te > Tjn- More recently and in this work,.-by Ipngei: sim- 
ulation time on larger systems, it is foundU-E£E3t£j that 
Te = Tra, and the VLL melts directly into the incoherent 
vortex liquid. 




FIG. ^ Top panel: Local Cooper-pair density < >, 
local condensate density | < ip' > p and helicity modulus 
along the field direction as functions of temperature for 
the Ginzburg-Landau model in a frozen gauge approximation 
with / = 1/60, F = 7, T'j,ip = 1, a^/^p = 6, and V = 60^ 
I < t/)' > P and Tz show a sharp drop at Tm, implying zero 
condensate density and lack of global phase coherence in the 
vortex liquid. 

Bottom panel: Same as top panel, but with T'j^jp = 0.3 
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ere have been sirmilations on the 



the Villain modeO, and the Lon- 
that suggest a two stage "melting transi- 
tion" . Here, the VLL melts at Tm into a coherent (disen- 
tangled) vortex liquid phase, where the free vortex lines 
are straight a well defined. In this phase, phase coher- 



B. Change in vortex-tangle connectivity 



We next discuss in some detail the results obtained 
for the quantity Ol, which probes the connectivity of 
the vortex-tangle in extreme type-II superconductors. 
We will make the following point: as for the zero field 
case, the increasingly sharp change in Ol from zero to 
one in finite field, with increasing system size, also de- 
notes a phase transition where a global [/(l)-symmetry 
is broken. This refers to a [/(l)-symmetry associated 
with the vortex-content of the Ginzburg-Landau theory. 
As argued in Section IIF, this symmetry of the vortex 
content of the theory is seen explicitly when rewriting 
it to a gauge-theory involving a local complex matter- 
field, see Section IIF. We also discuss the finite size ef- 
fects of Oli in systems with slab geometry, i.e. where 
Lx/Lz ~ Ly/Lz > 1, as well as in cubic systems. 

In Fig. ^ we show Ol as a function of temperature 
for several system sizes, V = 20'^, 60^, 120^. We see that 
for increasing system sizes, the largest temperature where 
Ol — increases, while the smallest temperature where 
Ol = 1 decreases. 

This finite-size effect in the behavior of Ol can not be 
explained by the 2D non-relativistic boson-analogy pic- 
ture. In this picture the transverse wandering u{z) 
vortex lines along the field direction (z )is given by''^ 



< {u{z)~u{0))'^ >=2Dz; D = 



(34) 



We see in Eq. ^ that the probability of finding a vor- 
tex line with finite line-tension traversing the system in 
a direction perpendicular to the magnetic field, without 
using PBC along the field direction, should decrease with 
increasing system size. Thus, if the vortex liquid regime 
is always describable as a liquid of vortex lines, then an 
inescapable consequence of this picture would be that Tl 
should shift to higher temperature with increasing sys- 
tem size. This is in clear contrast to the finite size effect 
of Ol shown in Fig. ^. In the zero field case, Ol = in- 
dicates that the line tension of vortex loops is finite, while 
= 1 indicates that the line tension of vortex loops is 
zero. The temperature where Ol jumps from zero to one 
is the critical temperature for the S-N phase transition. 
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FIG. 1^. Top panel: Ol as a function of temperature for 
the Ginzburg-Landau model in a frozen gauge approxima- 
tion with / = 1/60, r = 7, T'mp = 1, a^/^^ = 6, and 
V — 20^ , 60'^ , 120'^ . For increasing system size the largest tem- 
perature where 0_l = increases and the smallest tempera- 
ture where Ol = 1 decreases. Thus, in the thermodynamical 
limit, there exists a well defined temperature Tl where+Oi, 
jump sharply from zero to one. If we use the criteria Ol ~'C.9, 
to determine Tl, we find that Tl monotonically decreases to 
a limiting value for increasing system sizes. The inset shows 
the details of Ol close to Tl- Note how the curves for Ol all 
cross at the same temperature with increasing L. Note also 
how the lowest T at which Ol = 1 actually decreases with L. 
Bottom panel: Same as for top panel, but with T'j^ip = 0.3. 

We now focus on the inset of Fig. ^. Note how the 
curves for Ol cross, and reach, a value = 1 for pro- 
gressively lower temperatures as L increases. If a picture 
of the vortex-liquid in terms of well-defined vortex lines 
with non-zero linetension were applicable to this point, 
one would expect the point at which Ol reaches the 
value 1, to move monotonically up with L. The cross- 
ings of the curves for Ol observed in the inset of Fig. 



g simply would not occur. Note also the similarity of 
this finite-size effect, and the ones observed in Fig. [| 
for the zero-field case. There, it was argued that such a 
finite-size effect was strongly indicative of a percolation 
threshold for thermally induced unbound vortex-loops in 
the thermodynamic limitOu. The crossing point Tdoss 
seems a likely candidate for the limiting value of Tl as 
L ^ oo, see Fig. |l^ and the more detailed discussion 
below. This, in our view, provides strong numerical ev- 
idence that the progressively more abrupt change in the 
connectivity of the vortex-tangle as i — s- oo, is a real 
feature of the vortex system that survives in the thermo- 
dynamic limit, also at a finite magnetic field. In other 
words, the geometric transition signalled by the change 
in Ol seems to be a real feature and not an artifact of 
small systems. Whether or not it also corresponds to 
a finite-field thermodynamic phase-transition will be dis- 
cussed below. 

In the vortex representation, Eq. ^ the U{1)- 
symmetry to be broken is hidden, and can only be ex- 
plored implicitly using the conservation of N^. The con- 
nection is made explicit by rewriting the vortex Hamilto- 
nian in the disorder-field language, see EqjlJ of Section 
IIF. Below Tl, only field induced vortex line percolate 
the system. Thus, = Ny — and — fL^Ly. Here, 
fLxLy is the number of field induced vortex lines. For 
T > Tl, in addition to the field induced vortex lines, 
thermally excited "vortex lines" also exist. Thus, above 
Tl, Nf^ is not a conserved quantity and the global U{1) 
symmetry is broken, as for the zero field case. 

In Ref. 13 it was claimed that because the longitudi- 
nal superfluid density vanished precisely la^i-tiie melting 
line, as now found by several authoraiJiiJEj including 
the isotropic case, the vortex lines could not be consid- 
ered well defined in the vortex liquid phase. By itself, 
this is not a tenable conclusion. Nor does it follow au- 
tomatically that the vortex-lines are entangled arid, that 
the mechanism for VLL melting is entanglementca. To 
substantiate such a claim one has to investigate in more 
detail the geometric properties of the vortex tangle in the 
liquid phase, as done above and in Refs. 43 ij. Even if 
it should turn out that the loss of longitudinal superfluid 
density is entanglement it is probably more appropriate 
to view the entanglement as triggered by VLL melting 
transition rather than the converse. However, it is worth 
while pointing out at this stage that there is now con- 
sensus on the fact that at intermediate fields and above, 
the VLL melts into an incoherent vortex-liquid and that 
there does exist a regime where the molten phase consists 
of intact vortex lines, remarks to the contrary in Rcf. 
|T^ not withstanding. Moreover, various Monte-Carlo 
simulations agree that the Lindsnianaacriterion for VLL 
melting applies in this regimetSE^CJ'L^I. In the low-field 
regime, far less consensus has so far been reached. There- 
fore, the question of whether vortex loops influence VLL 
melting or not, and whether there exists a genuine transi- 
tion line Tl{B) inside the vortex liquid, are two separate 
issues. 
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FIG. |l^. Top panel: Helicity moduli Ta; and Tz, struc- 
ture function S{K.,kz = 0) and Ol for the 3-DXF-model as a 
function of temperature for / = 1/90, F = 7, and system size 
V = 72 X 80 X 80. 

Bottom panel: Tz and Ol for increasing system sizes. For 
increasing system size the largest temperature where Ol = 
increases and the smallest temperature where Ol = 1 de- 
creases. Thus, in the thermodynamical limit, there exists a 
well defined temperature Tl where Ol jumps sharply from 
zero to one, precisely as seen in the zero-field case. Shown is 
also specific heat for a system of size V = 360^ (shifted up by 
0.2A;_B for clarity). 

We hereafter focus on simulation results obtained for 
the iDXY-Tciodel. There is no quahtative difference 
between the results for this model, and the Ginzburg- 
Landau model. In Fig. |lO| we show results for the "iDXY 
model / = 1/90, F = 7. The top panel shows structure 
factor, superfluid density along the field, specific heat 
and Ol for a system of size 72 x 80 x 80. The bottom 
panel shows a sharpening of Ol for increasing system 
sizes. The trend in the change in the vortex-tangle con- 
nectivity is precisely the same as that seen for / = 1/60 
within the Ginzburg-Landau model including amplitude 
fluctuations. The lowest temperature at which Ol rises 
from zero, increases with system size, hut the highest tem- 
perature at which it reaches the value Ol — ^ decreases 
with system size. Again, we find a feature which indi- 



cates that a change in the vortex-tangle connectivity is 
undergoing a change. 



C. Effect of varying system aspect ratio 

According to the 2D non-relativistic boson-analogy of 
the vortex- liquid, Tl should be proportional to the aspect 
ratio L^/LjEi. 




0.9 0.95 1 1.05 1.1 



r 

FIG. Ol as a function of system temperature, obtained 
within the SDXY-model for / = 1/380, F = 7, for various 
system sizes. Top panel: Aspect ratio Ly/L^ — 1. 
Bottom panel: aspect ratio Ly/L^ — 2. Insets show details 
of the curve-crossings close to Tl- Note that while the lines- 
only approximation would predict a change in the crossing 
temperature of roughly a factor 2, they only change by about 
5%, which is within the uncertainty of the estimate for the 
crossing temperature. 

To further investigate the possibility of a breakdown 
of vortex-line physics inside the vortex-liquid regime, we 
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consider the crossing feature found in Ol in more detail 
for various aspect ratios L^/L^ of the systems on which 
the simulations are done. In this paragraph, we carry 
out the simulations using the 3DXY-mode\ with the pa- 
rameters / — 1/380 and F = 7. We have varied the field 
to illustrate that the features of are the same as for 
the higher fields, but do become sharper. Furthermore, 
comparing the results obtained from the FG model to 
the results obtained from the 3DXY-mode\, we again 
find that these models give qualitatively the same results 
when parameters are comparable. 
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FIG. ^ Finite-size effect in Tl as obtained for the 
3DXF-model using four sets of defining criteria for Tl, 
Ol = (0.10,0.20,0.90,0.95). Each of the criteria gives con- 
verging curves for Tl{L), whose limiting values are estimated 
by the crossing temperatures in the insets of Fig. |l^. Top 
panel: aspect ratio Ly/ — 1. 

Bottom panel: aspect ratio Ly/Lz = 2. Note that the hmit- 
ing values for Ly/ — 1 and 2 differ by about 5%, whereas 
according to a vortex-hne liquid picture, they should differ by 
about a factor of 2. 



In Fig. |TT|, we show Ol as a function of T for the 
3DXY-mode\ with parameters / = 1/380, F = 7, for a 
system of size Ly — L^ in one case, and L^ ~ Ly 

and Lz — 0.5Ly in the other case. Using the crossing tem- 
perature in the insets as an estimate for the temperature 
Tl in the thermodynamic limit, as for the zero-field case, 
we see that this temperature changes very little when 
changing the aspect ratio by a factor of 2. This indicates 
that in the thermodynamic limit there is only one Tl, 
regardless of the system shape. What these results in- 
dicate, is that the expectation one has based on the 2D 
non-relativistic bose-analogy of the vortex liquid, namely 
that Tl should scale with L^/L^, is not borne out. Note 
that the present case is very different from the situation 
encountered in the 3D Ising-model where a percolation 
threshold for overturned spins in an ordered spin state is 
found at a temperature which is lower than the critical 
temperatureL3. 

This may be further illustrated by considering the 
finite-size effect of T^, for two different aspect ratios 
Lx/Lz = 1 and Lx/Lz = 2. We investigate this by defin- 
ing Tl by four sets of criteria, namely the temperature 
at which Ol = (0.10,0.20,0.90,0.95). If the curves for 
Ol sharpen up, as seen in the above results, it is ulti- 
mately immaterial what sets of criteria are being used. 
The sets will give converging curves for Tl{L), one com- 
ing up from below and one coming down from above, see 
Fig. |l2|. We may use the best estimate for the cross- 
ing temperatures in Fig. |ll| as an estimate for what the 
limiting value of Tl will be in the thermodynamic limit. 
These results illustrate two important points, namely i) 
Tl does not move up monotonically with system size, but 
saturates at a specific value as L — > oo precisely as for the 
zero-field case, and ii) the limiting value of Tl is indepen- 
dent of aspect ratio. Both of these two points contradict 
expectations based on a vortex- Zzne liquid picture of the 
molten phase of the Abrikosov VLL. 



D. Scaling of the melting line T,n{B) 



In Fig. we show data from various simulations, 
of the vortex lattice melting line Tm{B). We want to 
emphasize the fact that there are two distinct scaling 
regimes for the melting line T,n{B), one at high fields 
which |We, somewhat arbitrarily denote high-field scaling 
regimdl23, and one at low magnetic fields which we iden- 
tify to be 3Z)XF-scaling. 

The dotted straight line is the curve given bytHl 
kBTm{B)/Jo = 0.41 y; < 4, where y = l/V/F. 
It describes the published numerically obtained melting 



lines for large enough filling fractions ?; < 4 or so well, 
in our case with F = 7 corresponding to approximately 
/ > 1/200. On the other hand, for y^ > 4, clear devia- 
tions from linear behavior is seen. 
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FIG. p^Top panel: Melting temperature T^{B) of the vor- 
tex lattice as a function of j/ = l/-y/7r. At large enough filling 
fractions, y < 2.0, Tl„{B) obtained from various simulations 
on the 3-DXy-model and boson-analogues of the vortex sys- 
tem, agree and is well described by T^{B) — 0.41y, the dotted 
line. At low filling fractions, y > 2.0, there is a crossover to 
3DXY critical-scaling of Tm(B)- The solid lines through the 
two sets of data points are 3DXY critical-scaling functions, 
described in the text. 

Bottom panel: Normalized melting temperature T!^{B) /T^ 
as a function of the varibale x/xq, where x — 1/%// and xq 
is an anisotropy dependent fitting parameter. Solid line is 
the 3DXY scaling function h(x) = 1 — [x / xq)^^^" , where 
u = 0.67. 

The melting curves obtained for F = 1 in Rcf . |l^ shown 
by the filled circles, and for F = 7 in Ref. |l^ shown by 
the half-filled circles saturate at low filling fractions / to 
the values given by the zero-field critical temperature, 
Tc. For F = 1, we have kcTc/Jo = 2.2, while for F = 7 
we have J^Tc/Jo = 1.12c3. The data given by the filled 
trianglealJ are obtained on the 3I?Xy-model with an 
anisotropy parameter F = 3. As F increases from 1, the 



zero-field transition temperature rapidly approaches 
its 2D value, although the transition is always 31? in 
character for finite anisotropy. Hence, the results from 
the anisotropic 3£'Ary-model with F = 3113 are very close 
to those of the 3DX Y-mo de\ with F = 7, see Ref. 13, 

The results of Ref. |l05| , obtained by fixing / — 1/15 
and var ying F G (1, ..10), agree entirely with our results 
of Refs.ll|,jl|,||in the low-field regime l/y/jT > 2. The 
significance of all these three sets of points is that they 
fall significantly below the straight line obtained from the 
high-field scaling of the melting line. 

Note also that even if we normalize the melting 
line T„i(BX quite arbitrarily, with a factor 1/(1 — 
Trn/7c)^''ll2J, this might take out the strong downward 
curvature of the data in the top panel of Fig. |l^, but 
there is absolutely no reason for why the slope of the 
resulting curve in the low field regime, which would be 
a straight line, should be the same as in the high-field 
regime. 

Assuming 3DXF-scaling for the melting line when 
y » 2, i.e. B/\l - T/T^l'^" = Bq where Bq is a field- 
scale that depends on anisotropy, we find kBTjn{B) / Jq = 
(kBTc/Jo) [1 — (xo/Fy)^/'^] on the melting line, and where 
the last term is negligible for low fields. Hence, we find 
that the melting line saturates to the true critical temper- 
ature Tc, as it obviously must. The dotted straight line 
kBTm/Jo = O.4I1;, overshoots Tc as the field is lowered. 
The Monte-Carlo results follow this line at large fields, 
but are however starting to be arrested in their tracks by 
the zero-field vortex-loop critical fluctuations already at 
around y = 2, thus crossing over to 3DXY critical scal- 
ing, as our Monte-Carlo simulations results show. 

In the top panel of Fig. n3, we have drawn the function 
Tm{B)/jQ - (T,/Jo)[l - (xo/yF)!/"] through the two 
sets of points obtained from Monte-Carlo simulations for 
y > A and F — 1,7, given by filled and half-filled circles, 
respectively. Using = 2.70 for F = 1 and xq = 6.45 
for F = 7, we find that the ZDXY scaling function given 
above fits the Monte-Carlo data well for y > A, while the 
high-field scaling is excellent for y < 2. Note how vastly 
different the scaling of Tm{B) in the two regimes y < 2 
and y > 4 is. 

The bottom panel of Fig. |l^ shows the low-field melt- 
ing line Tm{B) normalized by the zero-field critical tem- 
perature, obtained from simulations of the ZDXY model 
with F = 1, ksTc/Jo = 2.2 in Ref. P, F = 3, ksTc/Jo = 
1.34 in Ref. ^ and F = 7,kBTJJQ = 1-12 in this 
work and in Rcf. plotted in terms of the variable 
x/xo, where x ~ ^/Vf and xo is a fitting parame- 
ter for each F. The corresponding values of F and xq 
are (1,2.70), (3,4.65), and (7,6.45). For x/xq « 2 or 
less, i.e. at large enough fields, we see that deviations 
from Sil'X-scaling occur. For F = (1,3,7) this corre- 
sponds to 1// = (30,90,160), respectively. The line 
through the low-field data, is the 3DXY scaling function 
1 — {x/xo)~'^^'^ ■ Notice the sharp bending of the 3DXY- 
scaling function as x/xq increases beyond the value 3, and 
how the available numerically obtained melting curves 
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follow this line. This, in our view, provides strong nu- 
merical support for the notion that at low filling fractions 
/r^ << 1, the melting hne T„(B) obeys 3DXY critical 
scaling, while it follows follows a quite different "racan- 
field" type of scahng r,„(B) - 1/Vb at large ficldsH. 

The field above which deviations from 3DXY critical 
scaling is seen, increases as the anisotropy goes down. 
This is due to the fact that with increasing F, the melt- 
ing curve becomes flatter at low fields, leaving the critical 
region more rapidly as F increases. The width of the crit- 
ical region rappears to widen only slightly with increasing 
anisotropycS, whereas the melting line rapidly becomes 
flatter at low fields, such that this is the dominant effect 
in determining the field at which the melting line enters 
the critical region. 



E. Phase diagram, clean limit 

A summary of all of the above is contained in Fig.p^, 
we have included results from filling fractions 1/f G 
[90, ..1560]. The results we have obtained pertain to an 
extreme type-II superconductor in the absence of dis- 
order, since we are primarily interested in the intrin- 
sic properties of this phase-diagram excluding the se- 
vere complications due to disorder. There is a low- 
temperature vortex-line lattice phase. When the vor- 
tex lattice melts, it melts directly into an incoherent 
vortex-liquid with zero longitudinal superfluid density. 
The transverse superfluid density has been eliminated at 
temperature far below those where the VLL melts, by 
choosing low enough filling fractions to eliminate an un- 
wanted commensuration effect due to the presence of the 
numerical lattice on which the theory is defined. 

At zero magnetic field, we have demonstrated that an 
alternative way of describing the superconductor-normal 
metal transition, in addition to the phase-disordering pic- 
ture using the Ginzburg-Landau order parameter, is in 
terms of an unbinding of vortex-loops. We emphasize 
that although the quantity Ol we have focused on is not 
an order parameter, it may be tied to an order via the 
discussion in Section IIF. By including amplitude fluc- 
tuations explicitly in the Ginzburg-Landau theory, it is 
shown that this vortex-loop unbinding does not lead to 
critical amplitude fluctuations. A generalized "stiffness" 
characterizing the low-temperature phase which vanishes 
at the transition, is the long-wavelength vortcx-line ten- 
sion e(r), or equivalently the free energy per unit length 
of the the thermally induced vortex loops of the system. 

In a finite magnetic field, we find indications of a 
change in the vortex-tangle connectivity across the sys- 
tem at a temperature Tl{B), whose zero-field endpoint is 
Tc- This has been done by monitoring the quantity Ol in 
the Ginzburg-Landau theory or the 3DXY-m.odel in the 
same way as for the zero-field case. Ol has precisely the 
same characteristica at finite fields and zero field. In the 
regime Tm{B) <T < Tl{B), the connectivity across the 



system of the vortex-tangle of the molten phase is given 
entirely by the field induced vortex lines. This appears to 
change across the line Tl{B). We have been able to tie 
Ol to an order parameter even at finite field, see Section 
IIF, involving a breaking of a f7(l)-symmetry across the 
line Tl{B). 
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FIG. B — T phase-diagram for extreme type-II super- 
conductors based on Monte-Carlo simulations of the 3DXY- 
model with -B || c and F = 7. The phase-diagram splits into 
three differient regimes I, II, and III, characterized by the 
values of the Ginzburg-Landau and dual order parameters. 



At low magnetic fields, to the accuracy of our calcula- 
tions, we have found that the VLL melting line and the 
line Tl{B) merge at low fields. Below these low magnetic 
fields, the picture of the molten phase as a vortex-/me liq- 
uid appears to be questionable. For fields well above the 
point where Tm{B) and Tl{B) merge, we have found that 
the position of the VLL melting line is well described by 
a Lindemann criterion with Lindemann number c ~ 0.25, 
estimated from the Debye- Waller factor. 

Note that the rewriting of the theory Eq, 
is exact. The onset of the expectation value < 
place when vortex-loops unbind. Moreover 
Eq. Ilj exhibits an explicit ?7(l)-symmetry. When this 
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to Eq. |]J 

> takes 
the theory 



connection is made, it seems very reasonable to tie the 
observed change in the vortex-tangle connectivity to a 
vortex-loop unbinding and hence an onset of < >, i.e. 
the order parameter and the symmetry being broken in 
the transition, have been identified. 



VI. SUMMARY AND DISCUSSION 

We have explored the {B,T) phase diagram for ex- 
treme type II superconductors using two simplified ver- 



25 



sions of the the Ginzburg-Landau modeh i) The frozen 
gauge (FG) approximation where the gauge-field is fixed, 
while the phase and amplitude of the superconducting 
order parameter are allowed to fluctuate, and ii) the uni- 
formly frustrated 3D XY model where only phase fluc- 
tuations are allowed for. The former is obviously a more 
general model than the latter, while the latter is a com- 
monly accepted model in the studies of fluctuation effects 
in extreme type-II superconductors. Our results show 
that in the k oo-limit, where suppression of gauge- 
fluctuations is an exact feature of a superconductor, am- 
plitude fluctuations are completely dominated by phase- 
fluctuations over a sizeable temperature regime. The lo- 
cal order field < >, as well as the helicity modulus 
(global phase stiffness) develop an expectation value 
for T < Tc, and explicitly break the usual [/(l)-symmetry 
present in the Ginzburg-Landau theory. In contrast to 
this, the local Cooper-pair density < >, is finite 
both above and below T^. Our precise calculations close 
to T = Tc has brought out clearly its singular tempera- 
ture derivative at T = T^. Below, we list the main results 
of this paper. 
B ^ 

• In zero field, we have shown that the superconducting- 
normal metal phase-transition is described by a vortex- 
loop unbinding. This is achieved by correlating a de- 
tailed study of qualitative changes in the vortex-loop 
distribution function D{p) with calculations of super- 
fiuid density, condensate density, specific heat, ampli- 
tude fluctuations, and change in vortex-tangle connec- 
tivity, both including and excluding amplitude fluctu- 
ations of the Ginzburg-Landau order-parameter. The 
topological phase-fluctuations destroying the supercon- 
ducting phase-coherence are thus unambiguously identi- 
fied as thermally induced vortex loops. When amplitude 
fluctuations are included explicitly, they are found to be 
far from critical. In other words, the vortex-loop unbind- 
ing may not be viewed as a reparametrization of critical 
amplitude fluctuations of the Ginzburg-Landau order pa- 
rameter, as is sometimes claimed. 

• The vortex-content of the Ginzburg-Landau theory, 
formulated in SectionllE, is characterized by its own 
C/(l)-symmetry which becomes explicit on a further ex- 
act reformulation of the vortex sector in terms of a 
new gauge-field, see Section IIF. The low-temperature 
phase of the vortex-sector of the Ginzburg-Landau the- 
ory, where all vortex-loops are confined, exhibits a U{1)- 
symmetry. This symmetry of the vortex sector refiects 
the fact that there is a number conservation of vortex 
loops extending across the entire superconductor. In zero 
magnetic field, the conserved number is zero, and the 
distribution function for closed vortex loops of perime- 
ter p is an exponential function with length scale given 
by LoiT) = keT/eiT), where e{T) is the vortex-fine 
tension. At the zero-field critical temperature, we find 
vortex-loops with an algebraic distribution of perimeters, 
concomitant with an abrupt change in the connectivity 
of the vortex tangle in extreme type-II superconductors. 



The vortex line tension is found to vanish as a power law 
as Tc is approached from below, e(T) ~ |T — Tcp, with 
7= 1.45 ±0.05. 

• Both the change in the distribution function of closed 
vortex loops, and the abrupt change in the connectiv- 
ity of the vortex tangle, shows that there is a diverging 
length in the problem, i.e. Lo{T) — > oo; T —^ T~ . At 
this point, the number of closed vortex-loops extending 
through the system is no longer a conserved number equal 
to zero, it becomes finite and undergoes thermal fluc- 
tuations. Therefore, the C/(l)-symmetry characterizing 
the low-temperature vortex-phase is broken, the vortex- 
system has suffered a vortex- loop blowout, or unbinding. 

• The connection between the power-law behavior close 
to Tc of the vortex-line tension and the anomalous di- 
mension of the dual field (f) was discussed in Section 
IV C. Relating the power-law for the vortex- line tension 
to the susceptibility exponent 7 of the ^-ficld, in con- 
junction with the Fisher scaling law 7 = v{2 — rj^), al- 
lowed us to extract the value 77^ = —0.18 =F 0.07. This 
result was compared to renormalization group calcula- 
tion performed directly on the dual theory, for which the 
vortex-loop unbinding is the phase-transition, and excel- 
lent agreement was found. Note the negative sign of rj^ 
in the extreme type-II case. 

B 7^ 

• In finite field, we have studied the phase-diagram over 
a wide range of filling fractions /, corresponding to 
1// € [20, 1560]. Only a subset is shown explicitly, 
but all results are summarized in Fig. |lj. The VLL is 
found to melt in a first order phase transition, for all 
filling fractions considered, into a completely incoherent 
vortex liquid characterized by zero global phase coher- 
ence in all directions. At intermediate fields, the VLL 
melts into a liquid of vortex lines, whose position in the 
[B, T)-phase diagram is well estimated by the Lindemann 
criterion with a Lindemann number « 0.25. 

• We have performed a scaling analysis for the melt- 
ing line for all filling fractions considered. We find a 
crossover from mean-field type scaling at elevated fields 
to SDATy-scaling behavior at small fields, showing that 
for the aniosotropies we have considered, the melting line 
of the vortex-lattice at low fields is significantly affected 
by zero-field critical fluctuations in a sizeable region of 
the phase-diagram. 

• Significantly, in addition to the VLL melting tran- 
sition line Tm{B), we find indications of another transi- 
tion line, Tl{B), inside the vortex liquid. This line is the 
finite-field extension of the zero-field vortex-loop unbind- 
ing, and has an endpoint which is the zero-field critical 
temperature Tc- Below Tl, connectivity of the vortex- 
system is determined exclusively by the field-induced vor- 
tex lines. All vortex-lines threading the entire system 
are field-induced. Above Tl, this changes, as discussed 
in Section IVB. Above Tl{B) there exist vortex lines 
that thread the systen also perpendicular to the mag- 
netic field, without using periodic boundary conditions 
along the z-axis. 
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• We have performed a large-scale study of the finite- 
size effects in , and found that the temperature where 
the vortex-tangle connectivity changes does not move up 
with system size, like it would have done in a vortex- 
line liquid. In the 2D non-relativistic boson-analogy of 
vortex-liquid, such vortex-configurations are never found. 
The symmetry broken at Tl is a global U{1) symmetry, 
associated with the number-conservation of vortex paths 
threading the entire system, the considerations are simi- 
lar to the zero-field case. In a finite magnetic field, this 
symmetry is hidden in the usual Ginzburg- Landau local 
order field representation, but is brought out by a dual 
description of the Ginzburg-Landau theory. In zero field, 
Tl and are identical and there is only one phase tran- 
sition. 

• We have found that the vortex-system in the clean 
limit appears to be able to exhibit three distinct phases, 
I, II, and III shown in Fig. ^ characterized by the values 
of the Ginzburg-Landau order parameter < tp > and its 
dual order-parameter <(/)>. Here, we explicitly utilized 
the connection of Section HE between the vortex-tangle 
connectivity probe Ol and the t7(l)-ordering in the dual 
field (f>. We found the three regimes 

Region I : < tp >^ < (/) >= 0, 
Region U : < tp >= < cj) >= 0, 
Region III : < i/' < >^ 0. 

At low fields, we have found that region II vanishes. Note 
that the transition- line Tl{B) separating the regions II 
and III inside the vortex-liquid, was brought out solely 
through the dual description, it could not have been de- 
tected by studying the Ginzburg-Landau order parame- 
ter < ^ >, or any local function of it. 

A few further comments are in order. In the low- 
field regime, within a lines- only picture of the molten 
phase, one finds that the longitudinal correlation length 
of field-induced vortex lines above melting increases, due 
to the increased distance between field induced lines, be- 
ing given by 



(35) 



It was therefore pointed out in Ref. ^ that in order to 
correctly predict the direct transition from the Abrikosov 
vortex lattice to a phase-incoherent vortex liquid at low 
magnetic fields, or equivalently predict the direct transi- 
tion from the crystal phase to the superfluid phase of 2D 
non-relativistic bosons at T = at low magnetic, suf- 
ficiently large systems in the z-direction must be used. 
The use of too small systems could result in observing, 
merely as a result of a finite-size effect, a normal T = 
2D non-relativic bose-fiuid, or equivalently a disentan- 
gled vortex-liquid. The former cannot exist in the ther- 
modynamic limit in the absence of disorder, on quite gen- 
eral grounds. 

The above is a valid point of concern within the 2D 
boson-liquid analogy of the vortex system when looking 



for entanglement. It is no longer a point of concern if 
the lines-only approximation is abandoned and the con- 
nectivity of the vortex tangle is probed rather than en- 
tanglement. (Precisely how to establish a criterion for 
when field-induced vortex lines are entangled, also ap- 
pears to be problematic to say the least). For all fields 
we have considered, and for all sample geometries we 
have used, it is clear from our results that we have been 
able to correctly predict the direct transition from the 
Abrikosov vortex lattice to the phase-incoherent vortex 
liquid. The onset of and the change in the vortex- 
tangle connectivity is a separate matter. The vortex- 
configurations dominating the contribution to a change 
in Ol, are thermally induced unbound vortex loops and 
not field-induced fiux lines. Our results in the low-field 
regime are therefore not artifacts of considering too small 
systems in the z-direction. Quite the contrary, since we 
see the change in Ol also when making the system flatter, 
it supports the proposition that there exists a regime in 
the {B,T) phase-diagram, beyond the line Tl{B), where 
the notion of a vortex- Zme liquid physics most probably 
should be revised. 

The [/(l)-transition line Tl{B) has the zero-field su- 
perfluid normal state transition Tc as an endpoint. It is 
a feature of extreme type-II superconductors, even ho- 
mogenous, isotropic three-dimensional ones, and should 
moreover occur in Helium^ which is a perfectly three- 
dimensional, homogenous, isotropic superfluid. The pro- 
posed transition therefore is not in any obvious way 
connected to various previously proposed quite intrigu- 
ing scenarios leading to a loss of local line-tension of 
field induced— Ypiiices, often referred to as "decoupling- 



transitions" 
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. These phenomena rely on the layered- 



ness of the superconducting compounds, however they 
have no symmetry-breaking or order parameter associ- 
ated with them, but most importantly do not have a zero- 
field counterpart. Moreover, probing phase-coherence be- 
tween adjacent layers in aJavered superconductor as was 
for instance done in Ref.E£3 (see their Eq. 28) probes 
maximum g^-behavior, a part of reciprocal space not usu- 
ally associated with critical phenomena, which are in- 
frared singularities. Probing phase-coherence between in- 
creasingly distant layers ultimately amounts to comput- 
ing the helicity modulus Tj, which cannot vanish above 
the melting-line of the vortex lattice in the thermody- 
namic limit, in the absence of disorder. 

The Tl{B) line is potentially an important line in the 
{B,T) phase diagram. It locates the position in the 
{B, T)-diagram where the line-only approximation of the 
vortex-liquid breaks down. Pinning of vortices by ex- 
tended objects such as columnar pins may very well turn 
out to be inefficient beyond the line Tl{B). It also shows 
that the line-only approximation can be used to describe 
the vortex-liquid phase and the first order melting tran- 
sition of the VLL at Tm{B) only for large and inter- 
mediate magnetic induction. In low magnetic fields, on 
the other hand, Tl{B) and T,n{B) collapse into a single 
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linJi. Here, it would appear that a line-only approxi- 
mation does not describe the vortex liquid properly. The 
fields where the line only approximation fails in the en- 
tire liquidjegime are expected to be of order IT or less 
in YBColia. 
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